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ABSTRACT  OF  DISSERTATION 


NONLINEAR  METHODS  FOR  SPACECRAFT  ATTITUDE  MANEUVERS 

Spacecraft  attitude  control  is  described  by  a  nonlinear  dynamic 
mathematical  model.  When  Euler  parameters  are  used  to  define  the 
orientation  of  a  spacecraft  with  respect  to  an  inertial  frame,  the  model 
takes  the  form  ^  ^ s.  where  the  components  of  R.  are  quadratic  or 
higher  order  functions  of  the  state  variables.  Large-angle  maneuvers 
have  traditionally  been  accomplished  as  a  succession  of  small-angle 
rotations  from  a  sequence  of  defined  operating  points  about  which  the 
model  is  linearized.  Alternativately.  a  sequence  of  single  axis 
rotations  are  performed.  These  processes  are  computationally  burdensome 
and  are  unable  to  produce  fast  multi-axial  large-angle  rotational 
maneuvers . 

State-dependent  linearizing  transformations  allow  an  equivalent 
linear  system  of  the  form  x^  =  ^2  “  “*o-l  ~  *1-2  R  defined 
where  and  a^  are  freely  chosen.  A  control  strategy  u(t)  is 
determined  which  forces  the  linear  equivalent  system  to  evolve  in  a 
desired  manner.  The  input  £  to  the  physical  system  is  then  retrieved  by 
an  inverse  transformation.  Advantages  of  this  approach  to  spacecraft 
attitude  control  include  the  availability  of  linear  systems  techniques 
in  control  design,  realizability  of  fast  multi-axial  large-angle 
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maneuvers,  and  with  appropriate  choice  of  performance  criteria,  the 
ability  to  achieve  exact  terminal  states  by  specified  times. 

A  singularity  condition  exists  in  the  inverse  transformation  when 
detumbling  of  the  spacecraft  is  desired.  Sufficient  conditions  are 
developed  under  which  the  singularity  will  not  be  excited.  Imperfect 
knowledge  of  the  spacecraft  inertial  properties  leads  to  significant 
perturbations  from  the  terminal  state.  Introduction  of  a  linear 
regulator  acting  on  the  state  errors  in  the  linear  model  corrects  for 
these  perturbations.  Commanded  torques  generated  by  solution  of 
unconstrained  optimization  problems  for  the  linear  model  may  cause 
saturation  of  torque  actuators.  A  technique  of  reinitializing  the 
control  strategy  upon  encountering  saturation  is  introduced. 
Simulations  illustrate  the  feasibility  of  the  described  approaches  to 
spacecraft  attitude  control. 


Major  Albert  L.  Batten,  USAF 
Electrical  Engineering  Department 
Colorado  State  University 
Fort  Collins,  Colorado  80523 
Summer  1985 


iv 


ACKNOWLEDGMENTS 


A  successful  educational  program  cannot  be  achieved  by  an 
individual  without  significant  contributions  from  others.  I  wish  to 
acknowledge  some  of  the  people  who  have  made  mine  a  profitable 
experience.  The  United  States  Air  Force  Academy  sponsored  my  program, 
and  the  U.S.  Air  Force  Institute  of  Technology  made  it  financially 
possible.  Each  member  of  my  committee.  Dr.  Lee,  Dr.  Loeffler,  Dr.  Winn, 
and  Dr.  Zachmann,  have  had  a  direct  influence  on  my  work  outside,  as 
well  as  inside,  the  classroom.  Dr.  Dwyer,  my  adviser,  has  been  most 
helpful  through  the  entire  experience.  I  appreciate  his  patience.  Ms. 
Cheri  Emanuel  typed  (and  patiently  retyped)  this  manuscript.  A  person 
could  not  request  a  better  office-mate  than  Chen  Ning.  He  has  been  a 
constant  source  of  encouragement,  and  has  become  a  good  friend.  My  best 
friend,  my  wife  Nancy,  has  kept  our  family  together  through  this  period. 
To  say  I  love  her  does  not  begin  to  express  my  feelings. 

Above  all  else,  I  must  acknowledge  the  Lord  God,  who  created  the 
universe,  designed  the  natural  laws  that  it  should  obey,  and  gave 
mankind  wisdom  to  understand. 

Fort  Collins,  CO  A.  L.  Batten 
April  1985 


V 


TABLE  OF  CONTENTS 


Chapter  Page 

1.  INTRODUCTION .  1 

1.1  Background . . . 1 

1.2  The  Spacecraft  Attitude  Control  Model .  5 

1.3  A  Direct  Solution  Approach . 13 

2.  LINEARIZING  TRANSFORMATIONS .  19 

2.1  Bunt,  Su,  Meyer  Technique.. . 20 

2.2  Block  Triangular  Case .  29 

2.3  Freund  Technique .  34 

2.4  Summary .  41 

3.  OPTIMIZATION  PROBLEMS  AND  SIMULATIONS .  45 

3.1  The  Performance  Measure .  45 

3.2  Solutions  and  Simulations .  50 

3.3  Discussion .  61 

4.  DEALING  WITH  DIFFICULTIES  1 .  64 

4.1  Singularity  Avoidance . 65 

4.2  Variation  of  System  Parameters .  75 

5.  DEALING  WITH  DIFFICULTIES  II .  83 

5.1  Uncompensated  Saturation .  88 

5.2  Optimal  Decision  Strategy . . .  94 

5.3  Reinitialization  of  Optimization  Problem .  101 

5.4  Discussion . 106 

6.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY .  Ill 

6.1  Conclusions..... .  112 

6.2  Recommendations  for  Further  Study .  114 

REFERENCES .  117 

APPENDIX  A .  121 

APPENDIX  B .  125 

vi 


CHAPTER  ONE 


INTRODUCTION 


iji _ Background 

As  the  interest  in,  and  application  of,  space  vehicles  increase  in 
magnitude  and  complexity,  so  the  requirements  of  accuracy  and  timeliness 
of  control  strategies  for  those  space  vehicles  also  increase.  Our 
country's  present  administration  has  proposed  researching  the 
feasibility  of  a  space-based  defense  system  which  would  significantly 
reduce  our  vulnerability  to  long-range  missile  attack.  This  latest 
proposal  serves  to  highlight  the  advances  that  have  already  been  made  in 
space  technology  and  the  continuation  of  efforts  to  use  this  technology 
to  our  benefit. 

An  example  of  current  interest  is  the  Venus  Radar  Mapping  Mission 
planned  for  early  1986  [1].  The  objective  of  this  interplanetary 
mission  is  to  precisely  map  surface  characteristics  of  Venus.  The 
spacecraft  will  be  maneuvered  into  an  elliptical  orbit  about  the  planet. 
During  the  multiple  close  passes,  several  degrees  either  side  of 
pericention,  radar  images  of  the  surface  will  be  constructed.  As  the 
satellite  moves  toward  the  apocention  of  the  orbit,  it  will  be 
repositioned  to  transmit  these  images  to  receiving  stations  on  Earth. 
The  data  gathering  and  transmitting  sequences  will  require  extremely 
accurate  repositioning,  or  reorientation,  of  the  spacecraft. 
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Department  of  Defense  interest  in  high  performance  spacecraft 
attitude  control  is  evidenced  in  co-sponsorship  by  the  Air  Force  Office 
of  Scientific  Research^  with  NASA  Langley  Research  Center,  of  annual 
workshops  on  dynamics  and  control  of  large  space  structures.  One  of  the 
issues  emphasized  is  fast  maneuvering  of  these  spacecraft. 

"The  new  algorithm  research  has  been  driven  by 
requirements  for  spacecraft  control  software  that  is 
relatively  simple,  robust,  and  capable  of  being  executed  by 
prospective  on-board  computers  under  demanding  operational 
and  time  constraints...  Such  spacecraft  constitute  large, 
relatively  flexible  structures  that  must  slew  quickly 
through  large  angles,  then  rapidly  orient  and  stabilize 
their  attitudes  to  enable  their  on-board  acquisition, 
pointing,  and  tracking  (APT)  systems  to  acquire  their  next 
target .  "  [2] 


The  objective  of  current  research  in  spacecraft  attitude  control  is 
perhaps  best  summed  up  in  comments  made  by  Lt.  Gen.  James  A.  Abrahamson, 
Director  of  the  Space  Defense  Initiative,  to  Aviation  Week  and  Space 
Technology. 

"...;  what  we  have  to  do  is  move  our  own  system  very 
quickly  a  small  angular  distance  and  then  stabilize  and 
shoot  —  so  our  job  is  to  move  quickly,  stabilize  and  knock 
it  out  and  move  on  to  the  next."  [3] 


System  equations  for  multi-axial  reorientation  of  spacecraft  are 
nonlinear.  If  the  attitude  is  expressed  in  terms  of  Euler  angles,  then 
transcendental  functions  appear  throughout  the  model.  When  the 
spacecraft  orientation  is  defined  by  Euler  parameters,  the  system 
equations  are  functions  of  quadratic  terms  of  the  state  variables.  A 
standard  approach  to  dealing  with  this  nonlinear  model  is  to  define  an 
operating  point  about  which  the  model  is  linearized.  Small 
displacements  from  this  nominal  operating  point  can  be  effected  with 
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little  error.  The  magnitude  of  the  displacements  that  are  allowed  is 
governed  by  the  error  tolerance  acceptable  for  the  problem. 

For  a  multi-axial  large-angle  maneuver  to  be  completed,  successive 
operating  points  are  defined.  small  displacements  are  commanded,  and 
progress  along  a  prescribed  nominal  path  is  monitored  by  on-board 
sensors.  The  result  can  be  a  quite  accurate  reorientation,  but  this 
process  is  very  slow  and  computationally  burdensome. 

Single  axis  maneuvers  may  be  represented  by  linear  systems  models. 
An  alternative  to  the  linearization  process  described  is  a  sequence  of 
single  axis  maneuvers  where  only  one  reaction  wheel  or  thruster  pair  is 
used  for  each  maneuver.  A  time-consuming  reorientation  is  again  the 
resul t . 

To  circumvent  the  slow  successive  linearization  or  sequential 
single-axis  processes,  Junkins  and  Turner  [4]  formulated  a  nonlinear 
optimization  problem.  They  proposed  a  relaxation  process  for 
determining  costate  initial  values  which  then  allow  numerical  solution 
of  the  state  and  costate  differential  equations.  This  iterative  process 
must  be  performed  for  each  desired  set  of  initial  and  final  state 
values.  Solutions  determined  by  this  approach  converged  with  an 
accuracy  determined  by  the  specific  computer  in  use. 

Vadali  and  Junkins  [5]  suggested  using  the  method  of  particular 
solutions  to  solve  the  two  point  boundary  problem  for  state  and  costate 
histories.  As  presented  in  their  paper,  this  method  makes  use  of 
quasilinearization,  (linearization  about  a  nominal  trajectory  instead  of 
about  a  nominal  operating  point).  Again  an  iterative  process  is 
employed  to  converge  to  the  desired  solution,  and  storage  of  the  entire 
numerical  costate  trajectories  is  required. 


SECONDS 


SECONDS 


(a) 


Figure 


Attitude  Parameter  (y)  (b)  Angular  Velocity  (w) 


(c)  Applied  Torque  (t) 


1.3.  Simulation  of  single-axis  reorientation  using  only  linear 
feedback  terms. 
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3  2  2  3 

T  =  [6jY  +  ^2“^  '*’  ^5^“  ■*■  ^ 


The  first  two  coefficients  are  given  by 


62="  +  Q2*  Expressions  for  the  four  third  order 

coefficients  depend  on  the  closed  loop  transition  matrix  for  the  linear 
portion  of  the  state  equations.  They  are  lengthy  and  extremely  complex, 
however  they  are  achievable  with  the  use  of  simple  algebra  and  much 
determination  (see  App.  A). 

Simulations  were  run  for  I  =1.0  and  q,=q.=0.5,  the  same  values  used 

o  12 

in  [6].  State  trajectories  and  input  histories  (Figures  1.3  and  1.4) 
for  the  inclusion  of  linear  only,  and  linear  and  cubic  only,  terms  in 
appear  nearly  identical  to  those  in  Carrington's  example.  The 
advantage  of  the  Dwyer,  Sena  approach  is  the  closed-form  evaluation  of 
the  integral  expression  for  compared  to  the  solution  of  the  coupled 

differential  equations  in  the  Carrington  approach. 

Neither  method,  however,  is  without  its  difficulties.  The  amount 
of  computation  for  a  fairly  simple  spacecraft  maneuver  is  significant. 
When  multi-axis  control  is  desired,  computation  increases  substantially. 
Though  either  algorithm  shows  rapid  convergence  as  higher  order  terms 
are  included,  neither  approach  gives  an  exact  answer  to  the  optimization 
problem.  If  a  specified  (finite)  final  time  is  desired,  solution  of  the 
Ricatti  differential  equation  is  required,  which  is  still  feasible  in 
closed  form  by  the  method  of  Juang,  Turner  and  Chun  [21] .  The  higher 
order  terms  are  still  able  to  be  found  in  closed  form. 

A  more  appropriate  choice  of  orientation  parameter  for  the  single¬ 
axis  maneuver  is  the  angular  displacement,  6.  Under  the  assumptions  as 


outlined  for  this  problem,  the  model  would  be  linear 


u 


and  the  ^  are  n-homogeneons  nonlinear  terms. 


The  solution  can  then 


be  expressed  as 


T(t)  = 

where 


OD 

i  £.  +  X 

ss  2  •“  ss^  ^  ss 

ii=J 

"  ■  ^  *  « 

=  /  K^”^e«^i)dt 

o 

k=2 

2  ^^3  ss  ^  SS 

f  =  f  ~  y^p 

and  V  denotes  the  gradient  with  respect  to  For  the  problem  stated 

Af 

here,  f  is  a  constant  matrix  and  e^  is  easily  computed. 

The  nature  of  the  expansion  of  results  in  terms  appearing  in 

T(t)  that  are  of  odd  order  in  the  state  variables.  Hius.  when  only  the 

first  non-zero  term  in  J  is  used, 

ss 

t  =  +  62<*) 


and  when  the  second  non-zero  term  is  included 
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Although  present  technology  requires  external  thrusters  to  be 
operated  in  an  on/off  mode  only,  variable-thrust  gas  jets  are  currently 
under  development  [20].  Following  the  lead  of  Junkins,  Turner,  and 
Carrington  [4,6],  continuous  torque  maneuvers  are  designed  for  both 
internal  momentum  wheels  and  external  thrusters. 


1.3  A  Direct  Solution  Approach 

To  give  an  appreciation  of  the  amount  and  complexity  of  work 
involved  in  solving  the  nonlinear  problem  directly,  an  approach 
suggested  by  Dwyer  and  Sena  [18]  is  here  applied  to  the  single-axis 
model  (12a)  and  (12b).  Parameters  are  chosen  to  allow  for 
straightforward  comparison  to  a  similar  example  in  [6].  The  problem  can 


be  stated  in  a  format  compatible  with  Dwyer  and  Sena  as 


minimize  J  «  j  /  (£.^Q£.  +  T^)dt 
o 

subject  to  £  =  g.(Ji)  ^  V 

l(o)  =  (0.7071,  0.0)^ 

where  i.  =  (y»«i>)^»  =  (0,  Q  =  diag  (q^^.q^),  and 

o 


Expansion  of  £  (^)  results  in 

4  <t)  -  1 1  +  r 

n=3 


where 
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h(t)  =  i(r(t))  =  C(x(t))C(i.(tjj))“^h(t^) 
where  C{y_)  is  the  direction  cosine  matrix  written  in  terms  of  the  Euler 
parameters  and  the  vector  product  operation  x  such  that 

2  T 

C(x)v  =  (2y^-1)v  +  2yy  X  "  *  2^* 

Under  the  condition  ||x(t)||  <  1  for  all  t,  the  model  of  spacecraft 
attitude  control  with  internal  momentum  wheels  then  becomes 

X  =  (13a) 

u  =  (I  -I  )"^  (h(x)  X  i!i)  -  (I  -I  )~^  it.  (13b) 

Further  simplification  occurs  when  a  rest-to-rest  maneuver  is  to  be 

controlled  by  internal  torques  and  the  initial  wheel  angular  velocity  is 

zero.  Then  w( t  )  =  0,  fi( t  )  =  0;  therefore 
o  o 

4<‘„>  •  *  '.““o’  '  “ 

and  (13b)  reduces  to 

=  -  (I  -I  )'^  It. 
o  a 

Throughout  the  rest  of  this  study,  control  will  be  assumed  to  be 
either  external  thrusters  or  internal  momentum  wheels  but  not  both. 
Equations  (11a)  and  (11b)  or  (13a)  and  (13b)  are  the  respective 
mathematical  models  implied.  For  ease  of  reference,  a  general  nonlinear 
model  is  defined 


i  -  ii(t)  +  ¥(1)0. 

where  =  (i.j»£^)  =  (l.^»tt^)  •nd  £  is  either  i.  or  £  as  dictated  by  the 


selected  control. 
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This  reduced  aodel  is  included  to  permit  later  comparison  of  the  results 
of  Carrington  and  Jonkins  [6]  with  an  approach  by  Dwyer  and  Sena  [18]. 
Introduction  of  the  linearizing  transformations  will  also  begin  with 
this  model. 

If  the  spacecraft  orientation  is  to  be  controlled  by  the  momentum 
wheels  alone,  external  torques  are  zero  =  £)  and  total  angular 
momentum  in  an  inertial  frame  is  constant.  Junkins,  Vadali,  and 
Carrington  [S,6]  found  that  by  judicious  choice  of  another  reference 
frame,  this  assumption  allows  the  dynamics  equations  to  be  simplified 
and  the  order  of  the  model  to  be  reduced.  Dwyer  [19]  showed  that  given 
the  initial  value,  angular  momentum  in  body  coordinates  could  be 
propagated  as  a  function  of  the  Euler  parameters  by 
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i  =  r(x)tt 

*fexj4+T-]i  (8b) 

(I  ~I  )fi  =  ii  I  fe  -  I.  +  I  l“V  (8c) 

oa  oa 

Several  simplifications  can  be  made  for  specific  cases.  If  the 


momentum  wheels  are  locked  with  respect  to  the  spacecraft  and  control  is 
to  be  by  external  thrusters  only,  then  g  =  4  and  (8b) 

and  (8c)  become 

(I  -I  )w®Ijj)xia  +  T-  ji  (9) 

=  tt  X  -  1  +  it.  (10) 

Control  by  external  torques  alone  does  not  imply  ji  =  0.  Indeed,  from 
(10) 

k  =  X  tt  +  i). 

Inserting  this  expression  in  (9)  and  simplifying  gives 

lu^Iiaxw  +  T. 
o*  o 

Then  under  the  condition  ||x(i)||  ^  f  ^cr  all  t,  the  model  of  spacecraft 
attitude  control  with  external  jets  becomes 

X  *  nx)tt  (11a) 


(11b) 
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where  1^  is  the  matrix  of  momests  and  products  of  inertia  for  the 
complete  spacecraft  system  about  its  principal  axes,  is  a  diagonal 
matrix  of  moments  of  inertia  of  the  wheels  about  their  respective  axes 
(presuming  only  three  orthogonal  wheels  are  active),  and  fi  is  the  vector 
of  wheel  angular  velocities  with  respect  to  the  body  frame.  Inserting 
(3)  into  (2)  and  dropping  the  subscript  on 

Iji+lfi  +  ci)xh  =  T.  (4) 

o  a~  “ 

The  motor  torque,  ji,  required  to  change  wheel  angular  velocity  0  (and 
through  conservation  also  change  spacecraft  angular  velocity  gi)  is  given 
by 

H  =  I  (i  +  fi).  (5) 

a  “ 

Solving  for  I  Q  and  substituting  in  (4)  gives  the  equations  of  evolution 

a 

for  the  spacecraft  angular  velocities, 

(I  -I  )(fl  =  i  X  +  T  -  jl.  (6) 

o  a 

Similarly,  solving  (5)  for  jg  and  substituting  in  (4)  yields  the  wheel 
dynamics , 

(I  -I  )4  =  tt  X  -  I  +  I  I~^  11.  (7) 

The  development  of  the  dynamics  equations  given  here  parallels  that  in 
[171. 

The  system  model  is  then  given  by  (1)  for  the  kinematics  and  (6) 
and  (7)  for  the  dynamics,  repeated  here  for  completeness. 
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=  2  +  X  *  la) 

or  as  a  matrix 


Spacecraft  attitude  maneuvers  are  effected  by  the  application  of 
external  torques  (symmetric  pairs  of  jets,  or  thrusters)  or  internal 
torques  (momentum  transfer  wheels,  or  simply  momentum  wheels).  When 
three  wheels  are  used,  they  are  typically  aligned  so  their  axes  of 
rotation  coincide  with  the  principal  axes  of  the  spacecraft  (Fig.  1.1). 
If  a  fourth  wheel  is  included,  it  is  skewed  with  respect  to  the  other 
three  (Fig.  1.2). 

Dynamics  equations  for  a  rotational  maneuver  result  from 
application  of  the  law  of  conservation  of  angular  momentum.  The  rate  of 
change  of  angular  momentum  in  an  inertial  frame  is  equal  to  the  external 
torque  applied.  Because  torques  are  applied  along  the  spacecraft  axes, 
angular  momentum  is  referenced  to  the  body  frame,  but  the  rotation  of 
the  frame  must  then  be  accounted  for  [16], 

h.  (2) 
—inertial  “body  body  — 

where  h  is  the  angular  momentum  vector  referenced  as  shown  by  the 
subscripts,  and  x  represents  external  torques  supplied  by  thrusters. 
The  total  angular  momentxun  of  the  spacecraft  is  given  by 


(3) 


^3  ~  '*'  “2®^“  "'’  “3®^“ 

=  cos  ^(^)  +  (nj  +  n^  +  n^)siii^(^) 


Therefore  p^  is  related  to  £  by  the  expression 


Po  =  ±  \|1  -  1  . 

To  resolve  the  redundancy  inherent  in  the  full  Euler  parameter  set,  a 

»  9 

reduced  set  x.  i^  defined  by  x  =  fi  if  p^  >  0,  or  x  =  “  P.  if  p^  <  0,  and 
^  implied  constraint  in  this  definition  is  |  |x|  I  i. 

any  maneuver. 

As  the  spacecraft  orientation  changes,  so  do  the  angle  6  and  the 
unit  vector  n.  Whittaker  [IS]  shows  that  the  components  of  angular 
velocity  with  respect  to  the  body  frame  can  be  related  to  the  Euler 
parameters  and  their  time  derivatives.  Algebraic  manipulation  yields 
the  expression  for  the  time  rate  of  change  of  the  parameters. 


^o  =  2  <-^l“l  ■  ^2“2  - 
^1  =  2  ^^“l  ■  1'3“2  ^2“3^ 

^2  =  i  <l'3“l  ^  ^o“2  -  n“3^ 

^3=2  ^■^2“l  ^l“2  ^^0^3^ 


where  <j)  =  is  the  spacecraft  angular  velocity  vector.  The 

reduced  kinematic  model  may  then  be  written  as 


X  =  r(x)tt 


where  Rx)  can  be  considered  an  operator  such  that 


5 


in  soch  a  way  as  to  allow  specification  of  the  characteristics  of  the 
resulting  linear  equivalent  system.  The  application  of  interest  to 
Freund  was  3-  and  4-degree  of  freedom  robotic  manipulators. 

The  purpose  of  the  work  described  in  this  study  was  to  investigate 
application  of  the  theory  of  linearizing  transformations  to  the 
nonlinear  model  of  spacecraft  attitude  control.  When  a  linear 
equivalent  system  had  been  defined,  strategies  for  dealing  with 
variations  in  model  parameters  and  constraints  on  control  variables  were 
studied. 

The  remainder  of  Chapter  One  is  given  to  development  of  the 
spacecraft  attitude  control  model  and  one  direct  solution  approach. 
Three  algorithms  for  determining  a  linear  equivalent  system  are 
described  in  Chapter  Two.  Formulation  and  solution  of  optimization 
problems  for  the  linear  system  are  discussed,  and  simulations  are 
presented  in  Chapter  Three.  The  fourth  chapter  deals  with  parameter 
variations,  and  a  singularity  problem  that  arises  due  to  the  use  of  a 
reduced  set  of  Euler  parameters.  Compensation  for  commands  in  excess  of 
input  constraints  is  treated  in  Chapter  Five.  Recommendations  for 
further  study  are  presented  in  conclusion. 

1.2.  The  Spacecraft  Attitude  Control  Model 

An  arbitrary  orientation  of  a  spacecraft  with  respect  to  an 
inertial  frame  can  be  represented  as  a  single,  virtual  rotation  about  a 
unit  vector  n  which  is  fixed  relative  to  both  the  spacecraft,  or  body, 
frame  and  the  inertial  frame  [141.  This  representation  is  alternatively 

*r  »»r  •r 

described  by  a  set  of  Euler  parameters  £  =  (p  ,p. ,p.,p.)  =  (P  ,£  ) 

O  X  A  w  o 

where  p  =  cos  (^)  and  P.  ^  n. sin  (^),  i  =  1,2,3.  The  Euler 

O  ^  X  X  « 

parameters  are  not  an  independent  set,  for  from  the  definition. 


Carrington  and  Jnnkins  [6]  assumed  a  feedback  control  that  was  a 
polynomial,  rather  than  linear,  function  of  the  state.  This  approach 
requires  the  solution  of  a  matrix  Ricatti  equation  for  the  coefficients 
of  the  linear  feedback  terms.  Higher  order  coefficients  are  found  by 
solving  a  coupled  set  of  linear  nonhomogeneous  differential  equations 
which  depend  on  the  lower  order  coefficients.  Theory  was  presented  to 
include  either  fixed  or  free  final  time,  but  examples  only  illustrated 
the  steady-state  case.  The  number  of  terms  in  the  nonlinear  feedback 
required  for  satisfactory  convergence  is  apparently  a  function  of  the 
desired  maneuver.  For  the  examples  presented,  third  order  polynomials 
in  the  state  appeared  to  be  sufficient. 

Each  of  the  approaches  mentioned  solves  the  nonlinear  problem 
directly.  Each  depends  on  convergence  to  a  satisfactory  trajectory  of 
an  iterative  process  or  a  polynomial  series.  The  convergence  appears  to 
be  well-behaved  for  the  examples  presented,  but  the  result  must  be 
considered  to  be  approximate. 

Other  authors  have  sought  to  solve  nonlinear  control  problems 
indirectly,  by  in  some  manner  deriving  a  "linear  equivalent"  system. 
Improving  on  the  early  work  of  Brocket!  [7],  Hunt,  Su  and  Meyer  [8,9] 
proposed  the  derivation  of  a  nonlinear  change  of  variables  which  would 
transform  the  nonlinear  system  to  a  linear  system  in  Brunovsky  canonical 
form.  Sufficient  conditions  were  determined  for  assurance  of  the 
existence  of  such  a  transformation.  They  extended  the  theory  to  multi¬ 
input  systems  and  presented  examples  of  applications  to  flight  control 
systems  [10,11] . 

Freund  [12,13]  developed  an  algorithm  for  decoupling  nonlinear 
systems.  In  the  process  of  decoupling,  free  parameters  were  introduced 
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and  the  control  easily  determined.  The  point  of  presenting  the 
Carrington  and  Dwyer  approaches  in  the  context  of  the  single-axis 
maneuver  is  to  demonstrate  the  degree  of  difficulty  involved  in 
computing  the  solution  to  the  simplest  required  maneuver  when  the 
equations  of  motion  are  nonlinear. 

The  desire  then  is  to  find  a  way  of  formulating  the  multi-axis 
spacecraft  attitude  control  problem  that  gives  exact  results  for  either 
specified  or  unspecified  final  time  and  which  are  computationally  easier 
to  implement  than  those  described  above.  This  motivation  provides  the 
rationale  for  investigating  the  application  of  linearizing 
transformations . 


OIAPTER  TWO 


LINEARIZING  TRANSFORMATIONS 

The  objective  of  linearizing  transformation  theory  is  to  find  a 
change  of  variables  such  that  a  nonlinear  system  is  described  in  the 
transformed  space  as  an  equivalent  linear  system.  This  transformation 
will  be  nonlinear  itself,  and  the  equivalence  may  be  valid  over  the 
whole  or  only  a  portion  of  the  operating  universe.  A  generalized 
problem  can  be  stated 

given  £  =  (14) 


find  K  and  U 

such  that  X  =  X(£,),  u  =  !)(£., o) 

and  i  =  Ai  +  Ba  (15) 


for  some  state-  and  input-independent  matrices  A  and  B.  (For  the 
present  effort,  A  and  B  are  considered  constant.)  If  such  a  transform 
pair  X,n  can  be  constructed,  then  the  control  problem  can  be  formulated 
and  solved  for  the  linear  (£,n)  system  using  the  powerful  known  linear 
techniques.  The  input  £  is  then  transformed  to  nonlinear  (^,o)  system 
space,  and  the  maneuver  completed.  The  input  to  the  nonlinear 
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system  is  found  by  o  =  U  ^(i..u)  =  W(i,n)  where  the  inverse  exists. 
Usefulness  of  this  approach  is  eventually  determined  by  the  region  over 
which  D  ^  is  well  defined. 

As  in  linear  systems  theory,  control  can  be  implemented  in  either 
open-loop  or  closed-loop  form  as  'own  in  Fig.  2.1.  The  linear  problem 
can  be  solved  and  analytic  expressions  for  x(t)  and  u(t)  obtained.  The 
command  generator  supplies  not  only  the  linear  input  histories,  but  also 
the  linear  state  trajectories.  Both  are  fed  to  the  input  transformation 
W  to  produce  the  nonlinear  input  £  which  drives  the  system  (Fig. 
2.1(a}).  If  use  of  the  nonlinear  state  trajectory  in  determining  a  is 
preferred,  U  ^(^,u)  can  be  used  in  place  of  W(2.,n)  (Fig.  2.1(b)). 

Closed  loop  control  is  effected  by  transforming  the  state  variable 
to  linear  space  and  forming  the  linear  feedback  control  u  ==  u(x) .  The 
control  u  is  supplied  to  the  transformation  W  along  with  the  linear 
state  vector  (Fig.  2.1(c)).  As  before,  jj^  can  be  used  instead  of  x  by 
replacing  W  by  U 

Three  approaches  to  finding  an  equivalent  linear  system  will  be 
discussed . 

2.1  Hunt.  Su.  Mever  Technique 

In  their  paper  "Global  Transformations  of  Nonlinear  Systems"  [9], 
Hunt,  Su  and  Meyer  presented  sufficient  conditions  for  existence  of  a 
globally  equivalent  linear  system.  Their  technique  presumes  a  nonlinear 
system  of  the  form  (14)  with  a  scalar  input  o,  thus  the  matrix  V(j^) 
reduces  to  a  vector  field  j^(^) .  A  method  for  constructing  the 
transformation  between  nonlinear  and  equivalent  linear  systems  was 
introduced.  This  algorithm  is  dependent  upon  the  solution,  by  the 
method  of  characteristics,  of  a  system  of  linear  partial  differential 


(b)  Open  loop  control,  interface  driven  by  nonlinear  nodel  state. 


(c)  Closed  loop  control,  interface  driven  by  linear  model  state. 


Fignre  2.1.  Spacecraft  attitude  control  using  an  equivalent  linear 
system  for  input  determination. 
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eqnations.  When  the  sufficient  conditions  are  satisfied  and  the  system 
of  differential  equations  can  be  solved,  the  resulting  linear  system  for 
the  single  input  case  is  in  the  form  of  a  series  of  integrators. 


0  1  0  •  •  •  0 

0 

0  0  1  •  •  •  0 

0 

A  = 

• 

.  B  = 

• 

0  0  0  ...  0, 

1 

A  few  definitions  are  required. 

Let  ^(i.)  end  be  infinitely  differentiable  vector  fields  on  R°. 
The  Lie  bracket  of  with  respect  to  ^  is  also  a  vector  field  on  R°,  and 
is  defined  by 


Further,  define 


^  1  =  ^  gt  -  ji'  . 


(ad°  ^  t 

(ad^  io  =  u.  n 


(ad^  Hf)  ^  U,  (ad^"^ 

Finally,  let  be  a  function  on  R*^  which  is  infinitely  differentiable. 
Then  the  Lie  derivative  of  X^  with  respect  to  ^  is  defined 

dX.  dX.  dX. 

<dXi,ii>  =  ^  +  5^  ^2  ^  d^ 
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The  algorithm  for  a  single  input  is  summarized  as  follows; 

1.  Compute  (ad^  4,  ^  for  ‘  k=0, 1 , . . . ,n-l . 

2.  Construct  the  initial  curves  by  solving 

^  =  (ad“~^  4.  £.(0)  =  0 

to  get  £.(t^).  (n  equations) 

3.  Build  the  integral  surface  by  successively  solving 

^  =  (ad®"*^  jt.  f),  £.(tj,...,t^_^.0)  =  i.(tj....,t^_j) 

for  k=2,3....,n.  (n  equations  for  each  k) 

4.  The  solutions  to  the  partial  differential  equation  give  £.  as  a 

function  of  (t  j ,  t2. . . . .  t^^) .  Solve  for  tj^,  k=l,2,...,n  as  a 
function  of  • 

5.  Set  Xj=tj(£j^,£2»  •  •  •  solve  for  successive  X^'s  by 

\  =  <«iVi' 

6.  Find  U  by 

U  =  <dX^,  A  +  a>. 

Application  of  this  algorithm  to  the  single-axis  spacecraft  attitude 

control  problem  (12a  and  12b)  is  straightforward  and  clarifies  the 

T  T 

procedure.  Define  £,  =  ~  (t*w)  and  o  =  t.  Then  and  ^  are 

given  as 


Computing  the  necessary  Lie  derivatives. 
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(ad°  je) 


(ad^  A.  ^  =  U.  & 


The  system  of  partial  differential  equations 


1  I  2 

^  -  IT  0 

1  o’ 

5l(tj,0)  =  ;^j(tj) 

^2(0)  =  0 

Ii2(t^,0)  =  52(ti) 

is  solved  by  4j(tj,t2)  =  -  sinC^^^tj)  and  ”  ^^2* 

o  o 

easily  inverted  for  and  t2  as  ^1^1*52^  =  -  21 

*2^^1'^2^  =  Iq52'  setting  (i.)  “  =  -2I^sin  ^(5j). 

are  found, 

X2  =  <dX,.  <l>  =  -  1^52 


These  are 

sin 

X2  and  U 
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U  =  <dX2,  H  +  3{ib>  =  -o. 

The  negative  signs  that  appear  throughout  may  be  dropped  and  the 
resulting  transformation  in  terms  of  the  spacecraft  parameters  is 
defined 

X,  =  21  sin  y 

1  o 

X-  =  I  (I) 

2  o 

U  =  T. 


To  verify  that 
differentiate  x^ 


the  desired  linear  system  has  been  achieved, 
and  X2  vith  respect  to  time  using  the  chain  rule. 


X-  =  =  T  =  u. 

2  O 


Therefore 


i-[2 


u. 


As  mentioned  in  Chapter  One,  the  single  axis  control  problem  is  very 
quickly  solved  by  the  use  of  angular  displacement  and  angular  velocity 
for  the  state  variables.  The  linear  equivalent  system  was  developed 
here  to  introduce  the  use  of  the  Hunt,  Sn,  Meyer  technique  in  an  easily 
understood  context. 

A  later  publication  [101  extends  the  algorithm  outlined  above  to 
multi-input  nonlinear  systems.  The  objective  linear  equivalent  system 
becomes  a  Brunovsky  canonical  system  and  is  described  by  blocks  of 
series  of  integrators,  where  the  number  of  blocks  and  the  number  of 


inputs  are  equal.  Kronecker  indices  are  introduced,  specifying  the 
number  of  integrators  driven  by  each  input.  Multiple  inputs  add  little 
complexity  to  the  description  of  the  algorithm.  The  input  matrix 
is  partitioned 

Va)  =  [fjCf.)  I  I  •••  1 3^;a)] 

where  m  is  the  number  of  inputs.  Successive  Lie  brackets  of  HA]^)  with 
respect  to  calculated  for  i=l,2,...,m,  and  a  system  of  partial 

differential  equations  is  defined  as  before.  If  this  system  of 
equations  can  be  solved  for  ^  as  a  function  of  ( t^ , t2> • • • > t^) ,  and  if 
that  function  can  be  inverted,  then  the  transformation  is  constructed 
from  the  characteristic  parameters  associated  with  the  highest  order  Lie 
bracket  of  each 

Consider  now  the  model  presented  in  (11a)  and  (11b),  and  assume  a 
matrix  I  =  diag (I. , I. , I, ) .  Formatting  for  application  of  this 

O  .1  X  j 

algorithm  gives  i,  =  (i^.i^)^  =  Z  = 

0  0  0 

0  0  0 

0  0  0 

1/Ij  0  0 

0  l/I^  0 

0  0  l/I 

The  required  Lie  brackets  are  calculated  (see  Table  2.1),  and  the  system 
of  partial  differential  equations  defined. 
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Table  2.1.  Lie  Derivatives  for  the  Multi-axis  Control  Problem 


=  (0.  0.  0.  1/Ij.  0.  0)^ 


(ad^jt.jP^)  =  '^2  =  (0*  0*  0'  0* 


(adVf3)  =  ^3  =  <0-  0*  0'  0,  0.  I/I3) 


(-  ^  Ji  - 


.  _L  ,  0  bill  ,  ,  ,T 

21  M3'  21  M2'  M3'  1^X3  M2' 


(ad  ' 

M3’  “  HT  >  ^  I ^ 


_L. .  b:^ .  0 

21^  ^11'  Ijl2  M3'  I2I3  ^21' 


(ad  ^,5:^3)  =  [4.5^3]  = 


^  2I3  M2'  2I3  Ml' 


T  1,-1, 

|2  _3 _ 2 


"  21,  "  Ili-JI  '  I  I  M2'  1,1,  Ml' 
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.  £.(0)  =0 

U.  5^2^  •  =  i.(tj) 

[jt-  i^i3]  .  Ktj.t^.O)  = 

•  £.(tj»t2»t2,0)  =  i.(tj^,t2»t2) 

5^^  »  £.(  t j » t2 » ,  t^ . 0 )  =  £.( t j  ,  t2  *  tj  *  t^ ,  t ^ ) 

Because  of  the  stractore  of  ^(£.)  >  and  their  Lie  derivatives  for 

this  problem,  the  parameters  t^,  t2,  and  t^  will  appear  in  solutions  for 

and  not  in  solutions  for  i.2  •  Similarly,  t^,  t^  and  t^  will  appear 

in  solutions  for  £.2  alone.  Since  the  transformation  is  to  be 

constructed  from  the  parameters  associated  with  the  highest  order  Lie 

derivatives,  only  the  three  equations  involving  in  each  of  the  first 

three  sets  of  equations  need  to  be  solved.  (Note  that  in  general,  for  a 

2 

system  of  order  n,  a  total  of  n  partial  differential  equations  must  be 
solved  and  that  solution  set  inverted.)  If  the  nine  specified  equations 
can  be  solved  and  analytical  expressions  found  for  t^,  t2,  and  t^  in 
terms  of  algorithm  can  proceed  and  the  transformation 

constructed.  To  this  point,  that  problem  is  unsolved.  However,  the 
attempt  at  applying  the  Hunt,  Su,  Meyer  algorithm  to  spacecraft  attitude 
control  problems  is  not  unfruitful.  A  particular  strategy  for  choosing 


= 

dt. 


JLL 

dt. 


= 

dt. 


AL  _ 

dt. 


= 

dt. 


^  = 
dt. 
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the  linear  state  space  variables  is  suggested  by  the  integrator 
structure  of  the  desired  equivalent  linear  system. 


2.2  Block  Triangular  Case 

Combining  the  objectives  that  the  state  transformation  should  be 
input  independent  (X  =  the  input  transformation  should  depend 

explicitly  on  a  (U  =  and  the  linear  system  should  be  a  series 

of  integrators,  Dwyer  [19,22]  adapted  a  technique  of  Ileyer's  applicable 
to  nonlinear  systems  in  block  triangular  form.  The  technique  can  be 
nonrigorously  described  as  follows: 

Given  the  nonlinear  system  (14), 

1.  Find  a  set  of  elements  of  the  state  vector  ^  such  that 

(a)  the  total  of  the  selected  elements  and  as  many  time 
derivatives  as  are  independent  of  the  inputs  is  equal  to 
the  order  of  the  system,  and 

(b)  all  inputs  are  introduced  upon  the  next  time 
differentiation  of  the  selected  elements. 

2.  Define  I  to  be  the  set  of  elements  and  derivatives  found  in 
step  1(a). 

3.  Define  U  to  be  the  set  of  subsequent  derivatives  found  in  step 


The  single-axis  control  problem  again  serves  to  illustrate  the 
algorithm.  For  the  model 


=  7  \|l-y^  *» 


(I)  =  ^  T 
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n=2  and  two  elements  and/or  derivatives  must  be  selected.  Since  y  Y 
are  both  independent  of  the  input,  and 

••  djl  ■  .  ix  •  i.  2  1  .n  2 


explicitly  introduces  the  one  input  z,  that  choice  is  made.  Then  define 


Xj  =  Xj(y,«d)  =  Y 


*2  ~  X2(y>‘o)  =  Y  =  ^  \)  1-Y^ti) 


and  u  is  given  by 


n/  ^  12^  _1_  .r  2 

u  =  UlY.w.t)  =  Y  =  “  4  Yt^  +  Ml'Y 


The  resulting  linear  system  is  in  Brunovsky  canonical  form 


Also  important  to  note  is  that  a  requirement  to  drive  the  nonlinear 
system  to  the  origin  (YfW)  =  (0,0)  is  mirrored  in  the  linear  model  by  a 
requirement  to  drive  (xj,X2)  to  (0,0).  When  the  problem  is  solved  in 
terms  of  the  linear  system,  nonlinear  input  commands  are  obtained  by 
inverting  U(y.o),t), 


T  =  U  (y.w.u)  = 


(u  +  7  YU  )  . 


Knowing  the  desired  trajectory  for  x(t),  the  transformation  X  can  also 
be  inverted  to  find  desired  trajectories  for  y  u. 


T  =  *, 
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2  •  2 

<0  =  __  r  =  , - *2* 

A  singularity  occurs  in  the  expressions  for  t  and  u  when  y=1 .  But  the 
assumption  was  made  in  developing  the  model  in  Chapter  One  that 
I  Ixl  I  ^  single-axis  problem,  that  condition  translates  to 

-1  <  Y  ^  and  under  this  condition  the  singularity  will  not  be 

excited. 

This  method  of  finding  a  linear  equivalent  system  was  applied  to 
the  multi-axis  spacecraft  maneuver  controlled  by  internal  torques  [19], 
and  external  torques  [22].  In  both  cases,  choice  of  the  linear  state 
variables  follows  from  the  single-axis  problem  above,  and  the  desired 
result  is  achieved. 

For  the  external  thrusters  model  (11a)  and  (11b), 

-1  ^  =  X 

-2  ~  =  X  =  r(x)S!> 

u^  =  U(x.u.I.)  =  Y.  =  “  (  r(x)‘!>)X  +  ~  (  Rx)^)^ 

dx  du 

=  -  ■^(w^w)  +  Rx)I  ^  (I  wxti))  +  Rx)I 

(Linear  acceleration  associated  with  external  torques  will  be  designated 
Uj,  and  that  associated  with  internal  torques  will  be  designated  U2 
henceforth.)  The  state  transformation  can  be  inverted 

X  = 

w  =  r  ^  ^^^*2 


CHAPTER  THREE 


OPTIMIZATION  PROBLEMS  AND  SIMULATIONS 

3.1  The  Performance  Measure 

To  determine  a  desired  system  input  presuming  a  target  state 
normalized  to  zero,  the  modern  control  approach  defines  a  measure  of 
performance 


J  =  Jp(£.(t^))  +  /  [J^(£.(t),t)  +  J^(tt(t)  ,t)]dt 
^o 

which  is  then  minimized  subject  to  the  state  equations 

i  =  ^  (i)2. 

Typically  J  ,  J  and  J  are  chosen  to  be  quadratic  functions  of  the 
p  q  r 

state  and  input.  When  the  further  condition  of  linear  state  equations 
is  met,  solution  of  the  minimization  problem  reduces  to  the  solution  of 
a  matrix  Ricatti  equation.  In  many  cases  this  Ricatti  equation  can  be 
solved  and  an  optimum  control  strategy,  in  the  sense  of  minimizing  J,  is 
obta ined . 

However,  where  the  state  equations  are  nonlinear  the  complexity  of 
the  solution  is  greatly  increased.  The  nonlinear  nature  of  the  model 
for  spacecraft  attitude  control  (either  11(a)  and  11(b)  or  12(a)  and 
12(b))  results  in  the  inability  to  obtain  an  exact  solution  to  the 
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Free  parameters  in  the  generalized  transformations  allow  use  of  a 
combination  of  classical  and  modern  control  theory  to  be  applied.  The 
open  loop  linear  system  poles  can  be  placed  as  desired  by  appropriate 
choices  of  a^  and  a^.  Linear  inputs  n  can  be  determined  by  the  solution 
of  LQ  (linear  system,  quadratic  cost)  optimization  problems.  Chapter 
Three  presents  several  such  optimization  problems  and  results  of 
associated  simulations. 
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Table  2.3.  Generalized  Transformations  for  Spacecraft 
Attitude  Control  Using  Internal  Torques 


Nonlinear  Model 


X  = 

u  =  (I  -I  )~^(h(x)x«!))-(I  -I 
“  o  a  o  a 


=  i.  =  Rx)^ 

U2  =  T+a^i+a^X 


Inverse  Transformations 


X  = 

w  =  r~^(x)i2 


ii  =  -  2^(l!l^Ss!-4a^)  (I^-I^)x+h(x)x«+aj(I^-I^)!!)-(I^-I^)  r  ^(x)S2 


Linear  Model 


-1  “  ^2 

^2  =  ■  *0^1  “  *1-2  ^  *2 
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Table  2.2.  Generalized  Transformations  for  Spacecraft 
Attitude  Control  Using  External  Torques 


Nonlinear  Model 


i  =  r(x)l!> 

U  =  I  ^(I  (1>X(1))  +  I 
~  O  0~  ~  O  “ 


Transforma  t ion 


x^  =  i  =  r(x)<!> 

Uj  =  i+a^i+a^X  =  “  ^(a\)  +  r(x)lQ^<Io<!ixt!))+a2  nr)i!>+apX+nx)Io^l 


Inverse  Transformation 

X  =  ii 

“  =  r  ^ 

1  X 

T  =  z — (o)  <i)-4a  )I  Y-I  <i)xu-a,I  w+I 
~  2y^ - 0  0-^0 - 1  o- 


Linear  Model 


-1  -2 


-2  “o-r®l-2  -1 
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+  a,  r(i.)ia  +  a  Y  -  RaiXi  “I  )  ^ii.- 

1  —  o  o  a 

The  equivalent  linear  system  is  still  the  vector  form  of  a  damped 
second-order  system  vith  damping  coefficients  a^  and  a^  freely  chosen. 

2.4. _ Summary 

Although  the  algorithm  introduced  by  Hunt,  Su  and  Meyer  did  not 

give  an  equivalent  linear  system  for  the  multi-axis  spacecraft  attitude 

control  problem,  two  benefits  resulted  from  considering  their  work. 

First,  the  idea  of  even  looking  for  a  linear  equivalent  was  introduced. 

Second,  their  objective  of  a  linear  system  of  series  of  integrators 

suggested  a  starting  point  for  defining  the  transformations  that 

eventually  resulted.  The  two  difficulties  in  their  approach  are 

2 

solution  of  a  set  of  n  partial  differential  equations  and  the  inversion 
of  that  solution. 

Dwyer's  and  Freund's  techniques  differ  in  two  basic  ways.  Central 
in  the  Hunt-Su-Meyer  and  Dwyer  approaches  is  the  development  of  a 

transformation  that  will  give  an  equivalent  linear  system.  Freund's 
focus  is  on  finding  a  nonlinear  feedback  control  law  which  will  decouple 
the  system  outputs.  Secondly,  the  equivalent  linear  system  resulting 
from  the  first  two  approaches  is  predestined  to  be  a  series  of 
integrators,  while  the  linear  system  described  by  the  evolution  of 
Freund's  outputs  is  an  arbitrarily  damped  second-order  system. 

Application  of  Dwyer's  method  to  the  spacecraft  attituUe  control 
problem  can  be  easily  modified  to  achieve  the  generality  provided  b\ 
Freund.  In  place  of  u  =  x»  choose  u  to  be  n  =  x  +  a^X  instead. 

These  generalized  results  are  summarized  in  Tables  2.2  and  2.3  for  easy 


reference . 
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= 


.  V  = 


(I  -I  )  (h(x)»li))| 

L  o  a  J 


-(I  -I  )‘M 
L  o  a  J 


.  6  =  X. 


Calculation  of  N^(x)  N^(x)  >ce  the  same  as  before. 

9  9 

orders  d^,d2>d^  are  all  still  two,  although  for  this  case 


Differential 


d  ,.,1 


[N^X)]  ¥=  -  RxXi  -I  )‘ 

f  0  a 


Developing  Freund's  nonlinear  feedback 


8*  =  N^(x)  =  —  [Ni(x)]^  =  -  +  RxXI  -I  )  ^(h(x)iM) 

9  9  4  o  a 


A  =  -^  [Ni(x)]  V  =  -  r(x)(I  -I  )"^ 

9  o  a 


and  £.  is  unchanged.  Then  the  wheel  torque,  jj,  is  given  by 

iL  =  ^(x)  t~  4  +  RxXIp"!^)  ^(h(x)xw)  +  a^x  +  a^  P(x)so] 

-  (I  -I  )r"^(x>u, 
o  a  — 2 

which  simplifies  (somewhat)  to 


li  =  - 


1  /  T 


(w  «-4a  )(I  -I  )x  +  h(x)x<o  +  a,  (I  -I  Iw  -  (I  -I  )r  (x)]l, 
- o  oa  -  J-  —  loa  oa  — 2 


Definitions  of  X^(x>m)  and  X2(x>i!>)  are  unchanged  from  the  external 
torque  problem,  and 


£2  =  D(x.fi>.ii)  =  ■■  +  Rx)  ^(h(x)xii>) 
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IT  n  *  1 

=  - —  («  o)-4a  )I  X  “  I  ~  *,1  fi>  +  I  I  • 

2y^ - o  0  0  1  o  o  1 

Notice  specifically  that  a^=aj=0  gives  the  same  result  as  found  with  the 
Dwyer  approach.  Again,  the  transformation  is  defined  as 

=  XjCx.w)  =  X 

^2  =  =  i  =  HxJjii 

fij  =  =  X  +  a^x  +  a^jX 

=  “  ajRxJw  +  a^X  +  HxJIq^X* 

The  linear  equivalent  system  is  now 


■  0 

I 

O' 

X  = 

-a  I 
o 

1 

where  I  is  again  the  3x3  identity  matrix  and  will  remain  so  hereafter. 
As  before,  the  decoupled  nature  allows  system  analysis  and  controller 
design  to  be  completed  for  a  scalar  version,  then  generalized  to  the 
vector  case.  Evolution  of  x  =  Aj  ia  a  vector  equivalent  to  the  single¬ 
axis  result 

X  +  ajX  +  a^^x  Hi* 

When  Freund's  technique  is  applied  to  the  model  for  internal 
momentum  wheels,  much  of  the  work,  completed  for  external  thrusters  can 
be  reused.  The  elements  are  defined  by 
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N®(x)  =  X.  —  [N°(r)]  0 

^  di  * 

N^(;t)  =  [N°(r)]^  =  r(x)«!>.  lN^(x)l  ¥=  r(x)i~^  =*'  0,  ||xll<l 

*  dL  *  dL  ^ 


and  the  differential  orders  d^,d2«d^  are  all  two.  Construction  of  the 
elements  6  ;  A,  and  £  yields 


6*  =  N^(x)  =  ~  [Ni(x)]^  =  "  i  Xtjii^ia)  + 

"“fl  ^  oo 


A  =  -^  [Ni(x)]  r(x)l7^ 

dl  *  ° 


1 

£  =  X  *vN*(x)  =  a  X  +  *1  PCx^fit* 

kIC  0  o  * 

=o  ^ 


Here  the  assnnption  has  been  made  that  the  choice  of  free  parameters  is 

the  same  for  each  of  the  decoupled  channels.  If  that  were  not  the  case# 

a  and  a,  would  be  replaced  by  diagonal  matrices  A  and  A  ,  This 
o  1  o  1 

assumption  will  be  used  throughout. 

Nonlinear  input  commands  are  found  from 

t_=-A^[S.  +£]+A  ^u^ 

=  -  i^r"^(x)  I-  ^(w\)+r(x)i”^(ipS!>*“)+%x+02r(x)ii>]  +  ^(x)aj- 

”1  2 

Combining  terms  and  using  the  simplification  T  (x)x  = 

^o 

1  =  "  I  (z)e, 

2y  o  o  o  1  o  o  1 

O  V 
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Y  +  ajY  +  a^Y  =  Uj • 


The  equivalent  linear  system  is  found  by  setting 


*1  *  ^1  “  I' 


*2  ~  X2(y.w)  =  Y  =  i  \!  1-Y^« 


u,  =  D(y,w,t)  =  y  +  a,Y  +  a  Y 


1  2^1  .r  2  ^  ^  1  .f7" 

=  -  -  Yw  +2*1  “  *0“^  2I~ 


2 

Y  t 


so  that 


■  0 

1  ■ 

■O] 

X  = 

X  + 

.1  u, . 

-a 

L  0 

"*1. 

Ij  1 

Now  the  free 

parameters 

and  a^ 

can  be 

chosen  to  provide  desired 

second-order 

characteristics  to 

the  linear  system,  and  thus  to  the 

evolution  of  Y^t)*  A  special  case  is  a^^a^^^O.  Dnder  this  condition 
the  transformation  and  linear  equivalent  system  reduce  to  precisely 
those  found  using  the  Dwyer  technique. 

Setting  the  vector  output  6  =  x  multi-axis  problem  also  has 
the  desired  effect.  For  the  model  with  external  thrusters  (11a)  and 


(11b),  set  £.  =  (ij^,  ij)^  =  (X^.  as  before.  Then 


4(£)  = 


Rx)® 


r^(i^ttxa) 


0 

-1 


L^O 


Summarizing  in  vector  form  the  determination  of  differential  orders 
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the  components 


where  £.  =  (y.w 
chosen  8  =  y. 

nJ(i)  - 

Therefore  s 


and 


(d 


can  be  identified 


Ki)  = 


0 


2 

y  (i> 


^  il)  = 


Consistent  with  the  Dwyer  approach,  the  output  is 
First  the  differential  order  of  y  must  be  found: 

N°(y)  =  y.  [N°(y)]  4/ =  0 

^  dl  ^ 


1  \|l  -  y^m.  iNi(y)]  4/=  \|  1-y^  =>«  0  -l<y:<l. 

2  di 


;  2.  Construction  of  the  rest  of  the  elements  follows: 


8*  =  Nj(y)  =  -  iyu^ 

=  [Ni(y)]  ^  JlV  .  = 

dl 


21 


\[l^ 


p  =  r  =  “qY  +  *1(2 

k=o  ^ 


T  =  -  A  ^  [8*  +  p]  +  A  ^  u- 


21. 


I —  21 

[-  +  a^y  +  ai(2  yl~Y^“>J  ■*■  .  “j- 


Since  8  =  y  is  the  output,  the  algorithm  says  y  will  evolve  according  to 
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d.-l 

X 


*  i. 

£,(£.)  where  p.  =  X  a.  -N~  (6.)  is  the  i  element. 
1  .  i.k  4  1 

k=o  ^ 


Then  driving  the  nonlinear  system  with  o  given  by 


£=-A^t4  +fi.l+A^n 


results  in  decoupled  outputs  8^  that  evolve  according  to 


/j  ^  d.-l 
(d. )  1 

8.  ^  +  X  a.  .8.=u. 

1  . _  i.k  1  1 

k=o 


where  the  a^  ^  are  free  parameters  and  £  is  the  input  to  the  linear 
system 


^  =  A8  +  Bu. 


Notice  that  the  success  of  this  approach  depends  upon  the  existence 
of  A  Experience  (the  hard  task-master)  has  shown  that  the  wrong 
choice  of  outputs  can  result  in  a  A  matrix  which  is  not  invertible  at 
any  point  in  the  state  space.  No  discussion  is  given  by  Freund  to  guide 
the  correct  choice  of  outputs  when  that  freedom  exists.  Again  from 
experience  a  rule  of  thumb  seems  to  be  to  choose  a  set  of  elements  whose 
differential  orders  add  to  the  order  of  the  nonlinear  system.  This  is 
analogous  to  condition  1(a)  in  the  Dwyer  technique. 

For  one  last  time,  the  single-axis  maneuver  will  serve  to 
illustrate  the  procedure.  With  here  equal  to  the  scalar  moment  of 
inertia,  and 

•  i  Ti  2 
Y  =  2  “ 


2.3.  Freund  Tec 
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hnif^ne 

Freund's  [12,13]  objective  was  to  develop  an  algoritbm  to  yield 
nonlinear  control  laws  which  would  decouple  the  outputs  of  nonlinear 
systems.  A  side  benefit,  seemingly,  was  a  description  of  a  linear 
system  involving  the  outputs  and  their  derivatives.  The  description  of 
Freund's  technique  is  modified  to  highlight  the  result  of  interest  here, 
the  development  of  an  equivalent  linear  system.  In  fact,  this  algorithm 
as  presented  formalizes  the  non-rigorous  description  of  the  Dwyer 
technique  with  the  addition  of  freely  chosen  damping  terms. 

A  nonlinear  system  is  described 

i  =  |t(£.)  +  ^  (i.)2  ieR“,£eR“. 


Choose  m  elements  of  and  designate  them  6,  to  serve  as  outputs. 

Define  successive  Lie  derivative  operators  recursively  by 

(6.)  =  —  [N^“^  (6.)]^ 

^  "  di  ^ 

with 


(6.)  =6.. 

^  1  1 

Define  the  differential  order  of  6^  by 


d.  =  min  {j  :  —  ¥4  0,  j=l,2,., 

"  d^ 


.  ,n} 


The  algorithm  requires  the  construction  of 


•  *^i  th 

where  6.  =  N  (6.)  is  the  i  element, 
“  1^1 


d  ‘^•-1 

A(f)  where  A.  =  tN  ^  (6.)]  ¥  is  the  i 

"  di 


th 


row 
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Determinat ion  of  the  nonlinear  control  in  each  case  depends  upon  the 
existence  of  P  which  is  well  defined  everywhere  except  at 

=  0.  The  condition  =j*  0  is  equivalent  to  the  condition  |  |xl  I  ^  ^ 
which  has  been  introduced  before;  therefore,  when  the  operating 
assumption  is  satisfied,  U  exists  and  external  or  internal 

torques  for  the  nonlinear  system  can  be  determined. 

This  technique  results  in  identical  linear  equivalent  systems  for 
either  internal  or  external  control 

-[o 

T  T  T 

where  I  is  the  3x3  identity  matrix,  x  =  (x, ,x. )  ,  and  u  is  either  u,  or 

—  — 2  —  —1 

^2*  This  matrix  version  of  Brunovsky  canonical  form  is  equivalent  to  a 
set  of  three  decoupled  double  integrators.  The  decoupled  nature  allows 
controller  design  to  be  performed  for  a  scalar  model  and  the  result 
generalized  for  the  vector  case.  As  in  the  single-axis  scenario,  the 
origin  of  the  nonlinear  state  space  corresponds  to  the  origin  of  the 
linear  state  space.  Formulation  of  a  control  strategy  which  drives  the 
equivalent  linear  system  to  the  origin  will  have  the  same  effect  on  the 
nonlinear  system. 

One  drawback  tempers  (somewhat)  enthusiasm  about  the  result  of  this 
approach.  Poles  of  the  open  loop  linear  equivalent  system  are  at  the 
origin.  This  can  cause  a  stability  problem  as  will  be  shown  in  (Chapter 
Five.  The  third  technique  for  determining  a  linear  equivalent  system 
remedies  this  potential  problem. 
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and  as  before  the  nonlinear  input  can  be  expressed  either  as 
T  =  0  ^ or  as  t  =  WCi^.x^/U^).  The  former  is  notationally 
cleaner. 


T  =  ■7  I  r  ^  “  I  witu  +  I  r 

—  4o  o - o  - 1 


where  P  (x)  can  be  considered  an  operator  such  that 


r~^(x)Sj  =  ^  (YqJIi  +  -  YoY-iUj) 

or  as  a  matrix 

^o  V3  ^1^2  "^0^2  ^1^3 

-Y0Y3  +  Y1Y2  yl  +  Y2  Y^Yi  +  Y2Y3 

^0^2  ^1^3  "^0^1  ^2^3  ^o  ^  ^3  . 


A  useful  result  that  further  simplifies  the  expression  for  t  is 
—1  *y 

P  '*'(x)x  “  — X  •  Therefore, 

^o 

1  T  ”1 

T  =  r -  I  y(w  w)  -  I  U>XW  +  I  P  (x)]i,  • 

—  2r  o'^ -  o —  o  1 

o 

For  the  internal  momentum  wheel  model  (13a)  and  (13b),  the  state 
transformation  x  =  X(x»Sfi)  i®  left  unchanged.  However,  the  expression 
for  X  must  be  modified  because  of  the  difference  in 

U2  =  X  =  ~  (  P(X)“)X  +  ~  (  P(x)w)w 

dx  du> 

=  -  7  x(m\)  +  P(x)(I  -I  (h(x)x(o)  -  P(x)(I  -I 
4 -  oa  —  -^oa 

Nonlinear  control  a,  given  as  =  U  ^(x»S!i»fi)»  i® 
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optimization  problem.  This  difficulty  motivated  the  search  for  a  linear 
set  of  state  equations  that  are  equivalent,  under  some  change  of 
variables  to  the  nonlinear  model. 

If  an  optimal  control  is  now  sought  by  formulating  a  performance 
measure  which  is  quadratic  in  the  nonlinear  state  and  transforming  that 
measure  into  a  function  of  the  linear  state,  the  advantage  of  obtaining 
linear  state  equations  is  then  lost.  As  an  example,  consider  the 
problem  using  external  thrusters.  A  cost  function  is  defined 

*^f 

J  =  2  ^  2  +  q^w^C  t  )i!)(  t ) 

^o 

+  rT^(t)i(t))dt. 

Using  the  transformation  summarized  in  Table  2.2,  this  can  be  expressed 
in  terms  of  the  linear  state  variables  as 


t 


(qjX^(t)ij^(t) 


q2X2(t)r“^(ij(t))r  ^ij{t))i2(t) 


+  rt2^((oV4a^)I^3^  -  I^jjjxu)  -  ajl^u  +  I^T 


•  [2^^S!>^i!i“4ao)IoX  -  ^ 

*  O 
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where  r~^(x)  =  (T  No*  place  of  a  quadratic  performance 

measure  and  nonlinear  state  equations,  the  problem  has  become  one  of 
nonquadratic  performance  measure  and  linear  state  equations.  This 
hardly  represents  a  more  attractive  (or  easily  solvable)  problem. 

In  many  situations  the  objective  of  an  optimal  control  approach  is 
to  determine  a  control  strategy  that  will  drive  the  system  from  its 
initial  state  to  a  desired  final  state  with  bounded  inputs.  (Final  time 
may  be  either  finite  or  infinite,  i.e.,  specified  or  unspecified.)  The 
actual  form  of  the  cost  function  may  not  be  the  important  factor;  rather 
the  desire  is  to  find  a  solution  that  is  easily  computable.  When  this  is 
the  case,  the  performance  measure  can  be  formulated  as  a  quadratic 
function  in  the  linear  state, 

^f 

J  =  ~/(tj)Pji(t^)  +  ^/  (i’'(t)Qi(t)  +  u^(t)Ra(t))dt,  (16) 

^o 

and  then  minimized  subject  to  the  linear  state  equations.  The 
simplicity  of  the  equivalent  linear  system  developed  in  the  previous 
chapter  yields  several  situations  in  which  an  exact  analytic  expression 
for  the  "optimal"  control  can  be  found. 

For  the  spacecraft  attitude  control  problem,  there  are  two 
important  conditions  which  most  be  satisfied  so  that  a  solution  to  the 
linear  model  optimization  problem  will  transform  to  an  acceptable 
solution  to  the  nonlinear  problem.  First,  a  bounded  solution  for  u 
transforms  to  a  bounded  solution  for  o  only  if  Yq  ^  0  (i.e., 


This  can  be  accounted  for  in  the  problem  statement  by  the  introduction 


However,  this  nonlirear 


2  2  2 

of  an  inequality  constraint  ■*’*12^*13^^* 
constraint  introduces  difficulty  in  solving  the  problem  just  as  the 
nonlinear  state  equations  did.  Another  approach  is  to  solve  the 
unconstrained  problem  in  general,  then  determine  a  set  of  initial 
conditions  for  which  the  solution  results  in  bounded  torques.  I>is 
approach  is  discussed  further  in  the  Chapter  Five. 

A  second  requirement  is  that  the  terminal  state  in  the  linear  model 
must  transform  to  a  unique  state  in  the  nonlinear  model.  There  is  one 
state  which  is  preserved  through  both  the  developed  transformation  and 
its  inverse.  The  only  state  in  the  nonlinear  model  that  is  associated 
by  the  transformation  with  the  origin  in  the  linear  model  is  the  origin. 
If  the  terminal  state  for  reorientation  maneuvers  can  always  be  defined 
as  the  origin,  there  will  be  no  ambiguity.  Driving  x.  to  the  origin  is 
equivalent  to  aligning  the  principal  axes  of  the  spacecraft  with  the 

inertial  axes.  Since  the  inertial  frame  can  be  arbitrarily  chosen,  the 

target  state  can  always  be  the  origin.  Further,  ^2=0  implies 
w  =  r  ^ (x)i2  *  — *  Therefore,  when  detumbling  is  desired  (final 
orientation  arbitrary),  a  solution  which  drives  x^  to  0  will  result  in 
the  spacecraft  coming  to  rest. 

There  is  yet  some  flexibility  in  constructing  the  performance 

measure  (16).  Final  time  can  be  either  finite  or  infinite.  Choice  of 
finite  final  time  results  in  a  trajectory  which  will  be  within  a 
neighborhood  of  the  desired  terminal  state  by  time  t^.  The  size  of  the 
target  area  is  specified  by  choice  of  .  Exact  terminal  control  is 
obtained  when  ®  (that  is,  the  diagonal  elements  of  P^  become 

infinite)  [23].  Infinite  final  time  results  in  a  trajectory  which 
asymptotically  approaches  the  desired  final  state.  Although  the  target 
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state  need  never  be  reached  in  finite  time,  there  is  one  reason  why  the 
final  time  might  be  made  infinite.  Closed  loop  control  with  t^  finite 
results  in  time  varying  feedback  gains.  When  t^  — >  the  steady  state 

matrix  Ricatti  equation  is  solved,  rather  than  the  time  varying 
equation.  Resulting  feedback  gains  are  then  constant,  and  closed  loop 
control,  with  the  associated  benefits  of  state  feedback,  can  be  more 
easily  used. 

Accepting  the  premise  that  the  objective  of  the  optimal  control 
problem  formulation  is  to  design  a  trajectory  that  can  be  easily 
computed,  the  state  penalty  matrix  Q  can  be  set  to  zero  when  final  time 
is  finite.  This  allows  solution  of  the  time  varying  Ricatti  equation 
for  a  gain  matrix  P(t)  to  be  replaced  by  solution  of  a  time  varying 
Lyapunov  equation  for  P  ^(t).  Only  first  order  terms  appear  in  the 
Lyapunov  formulation  and  the  solution  is  correspondingly  easiei.,  If  Q  ^ 
0,  a  solution  to  the  time  varying  Ricatti  equation  may  be  obtained  using 
a  method  (somewhat  more  involved)  suggested  by  Juang,  Turner  and  Chun 
[21]. 

When  the  final  time  is  infinite,  a  positive  semi-definite  Q  will 
force  the  state  to  asymptotically  approach  the  desired  terminal  state. 
Then  P^  may  be  set  to  zero  and  the  steady  state  Ricatti  equation,  which 
is  quadratic  in  the  matrix  P,  must  be  solved.  No  benefit  is  derived  in 
this  case  from  trying  to  solve  for  P  ^  rather  than  P  for  the  non-zero  Q 
will  make  the  equation  in  P  ^  quadratic  as  well.  With  P^=0,  the  state 
penalty  matrix  must  not  also  be  set  to  zero  for  then  the  performance 

OD 

1  T 

measure  would  become  J  =  ~  /  u  jidt.  With  no  penalty  on  the  state 

o 

trajectory  or  the  terminal  state,  the  optimal  solution  is  u=Q^  and  the 
system  will  remain  in  its  initial  state. 
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Some  advantages  can  be  realized  by  including  a  penalty  on  u  [24]. 
This  is  accomplished  by  defining  an  augmented  state  vector  x  = 

(1)  • 

increasing  the  order  of  and  Q  appropriately  and  replacing  u  by  u  =u 

in  (16).  In  this  situation,  the  equation  u=u^^^  is  added  to  the  state 
equations  and  initial  and  final  values  for  u  may  be  freely  chosen. 

3.2.  Solutions  and  Simulations 

With  the  equivalent  linear  system  developed  in  the  previous  chapter 
and  the  foregoing  discussion  of  performance  measures  in  this  chapter,  an 
infinite  number  of  linear-quadratic  (LQ)  optimization  problems  can  be 
formed.  The  damping  terms  must  be  picked;  choice  of  specified  or 
unspecified  final  time  must  be  made  and  if  specified,  the  value  must  be 
decided.  For  finite  time  problems  must  be  chosen;  for  infinite  time 
problems  Q  must  be  chosen;  and  the  choice  of  whether  or  not  to  penalize 
u  must  be  made.  For  the  purposes  of  example  simulations,  a  desired  type 
of  maneuver  must  be  defined  with  appropriate  initial  conditions  and 
terminal  states.  Control  actuators  (wheels  or  thrusters)  must  be 
chosen,  control  strategy  (open  loop  or  closed  loop)  must  be  picked,  and 

the  spacecraft  parameters  (I  and,  if  wheels  are  used,  I  )  must  be 

o  a 

specified.  The  most  general  formulation  of  the  problem  statement  is 
given  by 

^f 

f  (q|X^(t)Xj (t)  +  q2Z2(t)i2^i) 


1 

^  2 
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+  +  q^u^(t)u(t))dt  (17) 

subject  to  ij  =  *2  ~  ~1  — ”  “1 

£-2  -  -2^^^  “  -2  “  -2’ 

From  this  universe  of  possibilities,  five  specific  cases  are 

selected  to  illustrate  the  use  of  linearizing  transformations  in 

spacecraft  attitude  control.  The  maneuver  is  selected  to  allow 

comparison  to  the  results  of  Vadali  and  Junkins  [5].  A  slowly  tumbling 

spacecraft  is  to  be  stabilized  at  a  specified  final  orientation. 

T 

Initial  conditions  are  =  (0.442,0.442,0,442)  and  (^(0)  = 

(0.01,  O.OOS,  0.001)  with  the  terminal  state  being  the  origin  as 
discussed.  The  inertial  properties  are  described  by 
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0.00 
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off-diagonal 

terms  in 

the  system 

inertia  matrix  are 

due 

to  the 

presence  of  a  spare  skew  wheel  which  is  not  active. 

For  the  purpose  of  these  illustrations,  unconstrained  input  torques 
are  supplied  by  three  momentum  wheels  with  axes  along  the  principal  body 
axes.  Perfect  knowledge  of  the  system  parameters,  dynamics,  and  state 
and  no  external  disturbances  are  assumed.  This  allows  an  open  loop 
control  strategy  to  be  used  in  each  case.  However,  in  the  solutions 
presented  the  closed  loop  control  law  is  always  found  first  and  could  be 
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used,  when  appropriate.  When  finite  final  time  is  specified,  — >  »>  so 

that  exact  terminal  control  is  achieved.  The  remainder  of  the 

specifications  are  identified  as  each  case  is  discussed. 

Whether  damping  coefficients  are  chosen  identically  for  all  three 

channels  so  that  ~  *^1—2  *  chosen  differently  for  each 

channel  so  that  x  =  -  A  x,  -  A.i.  +  u  where  A  and  A,  are  diagonal 
~  o“l  1“2  —  o  1 

matrices,  the  equivalent  linear  system  is  decoupled.  Then  the  entire 
equivalent  linear  system  can  be  treated  as  the  combination  of  three 
scalar  second-order  systems.  Accordingly,  optimization  problems  are 
solved  in  the  scalar  case  and  generalized  to  the  vector  system. 
Identical  desired  response  characteristics  are  assumed  for  each  channel. 
When  the  damped  system  is  employed,  a^  and  a^  are  chosen  to  provide  a 
critically  damped  response  and  are  the  same  for  each  scalar  system.  A 
summary  of  the  five  examples  follows  with  specifications  referenced  to 
(17).  Detailed  solutions  to  the  five  problems  presented  are  found  in 
Appendix  II. 

3.2.1  Case  I 

Finite  final  time  is  selected  as  t^=100  sec.  The  norm  of  the 
derivative  of  the  input  u  is  not  penalized,  therefore  q^=0  which  implies 
p^  must  also  be  zero.  Terminal  penalties  p^  and  p^  are  infinite.  No 
penalty  is  assumed  on  the  norm  of  the  state  throughout  the  trajectory. 


so  ^arefore,  the  performance  measure  is  given  by 


53 


and  is  allowed  to  become  infinite.  The  system  equations  describe  a 
double  integrator,  therefore  a^=aj^=0.  Solution  of  this  optimization 
problem  yields  a  closed  loop  control  law 

u^Ct)  =  -6(t^-t)~^  ij(t)  -  4(tj-t) 

Analytic  solutions  can  be  found  for  z^(t),  Z2(t)  and  u(t) 

ij(t)  =  c.j(t£-t)^  + 

i2(t)  =  -  3cj(t^-t)^  -  2c2(tj-t) 

=  6£j(t£-t)  +  2£2 

where  =  St^  These  expressions 

are  processed  through  the  inverse  transformation  for  use  with  momentum 
wheels,  and  resulting  histories  for  y^,  u,  and  £  are  shown  in  Fig. 
3.1(a)  through  (d). 

3.2.2  Case  II 

With  t^  still  100  secs  and  a^  and  a^  remaining  zero,  a  penalty  is 
now  placed  on  | |u| | .  Therefore,  Pj,  P2.  and  Pj  are  all  infinite  and 
<lj^,  q2»  and  are  all  zero.  The  performance  measure  is  now  given  by 

J2  =  f  x^(tj)P^i(tj)  +  ^  /  (i^(t)Ri(t))dt 

o 

where  again  P^  — >  «d.  This  problem  results  in 

i2(t)  =  -60(tf-t)"^£j(t)  -  36(tf-t)“\2(t^  ~  9(t£-t)“^u(t) 

which  must  then  be  integrated  before  being  fed  through  the  inverse 
transformation  and  applied  to  the  nonlinear  system.  Expressions  for  the 
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(c)  Angnltr  Velocity  (jg) 


(d)  Applied  Torque  (^) 


Figure  3.1.  Stabilization  and  reorientation  minimizing  performance 
measur'^  using  momentum  transfer  wheels. 

D  >  1st  component,  v  >=  2nd  component,  0«  3rd  component 
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the  evolution  of  x^,  X2  and  u  are 


ij(t) 

i2(t) 

u^d) 


=  £2(1^-!)^  +  £2(t£-t)  +  £2(t£-t) 

=  20c2(t^-t)^  +  12c2(t^-t)^  +  dCgCt^-t) 


1.  o  .  ,.-4  o  n  «*“3  o  ,-^~4  o  o  -2  o 

where  =  6tj.  x^  +  3t^  x^  +  0.5t^  u  ,  £2  “  “ISt^  x^  -  7t^  *2  ~  ' 

and  c_^  =  lOt^^Xj  +  4t^^X2  ■*"  0.5t2.^n°.  Simulations  for  the  nonlinear 
variables  are  in  Fig.  3.2(a)  through  (d). 


3.2.3  Case  III 

The  third  example  is  an  infinite  time  case.  As  discussed  earlier. 
P2<  and  p^  are  set  to  zero.  Penalizing  the  norm  of  u  implies 
and  q^=0.  State  penalties  and  42  are  nonzero  to  drive  (i2,X2)  to  the 
origin.  They  are  picked  by  trial  and  error  to  be  =  9x10  ^  and 

_3 

42  =  6x10  to  make  the  x  trajectory  approximate  that  found  in  case  I 
and  case  II.  Performance  measure  is  described  by 

CD 

J3  =  2  /  (i^(t)C(s(t)  +  u^(t)u(t))dt. 
o 

Damping  terms  a^  and  a^  are  left  zero.  The  feedback  control  is  given  by 

%  -  -  \Fr  *1  -  i2- 

State  trajectories  and  the  open  loop  control  law  are  represented  as 


/. \  .  ~“t  ,  -ot 

Xj(t)  =  £2te  +  £2e 

X2(t)  =  -a£2te  +  (£2  -  <*£2)® 

u^Ct)  =  o^£2te  -  (2o£2  -  a^£2)e 


WI«N5  /  SCCONO 
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where 


£j  =  ax°  ^2'  ~2  ~  ~l'  ~  •  Figures  3.3(a)  through 


(d)  give  the  resulting  curves  for  the  nonlinear  system. 


3.2.4  Case  IV 

An  infinite  time  problem  with  a  penalty  on  |  |u|  |  is  solved  in  case 
IV'.  The  performance  measure  in  this  case  is 


=  ^  /  (i^(t)Qx(t)  +  u^(t)u(t ) )dt . 
o 


Terminal  penalties  and  damping  coefficients  are  all  zero.  State 

penalties  <l2  acceleration  penalty  are  again  chosen  by 

-7  -4  -2 

trial  and  error,  q^  =  1x10  ,  q^  =  3.23x10  and  q^  =  1.39x10  ,  and  q^ 

is  set  to  1.  The  closed  loop  expression  for  u  is 


n^Ct)  *  -  -  P23i2^^^  “  ^33*3^^^ 


where 


=  \pr- 


satisfies 


-  2^2^23 


**^1^23 


+  ^2  ■  ‘*‘*l‘*3^  =  ®  P33  = 

q^  are  picked  here  such  that  p^ 


|2p22+q2 .  (The  values  for  q^,  q^,  and 
P —  3  2 

=  Jqj  =  a  ,  p^j  =  3a  ,  and  p^^  =  3o.) 


When  u  is  integrated  and  applied  to  the  linear  system,  the  results  are 


Xj  (t ) 

X2(t) 

U2(t) 


2  -at  .  .  -ot  .  -at 

c.jt  e  +  £2^®  — 3® 

-oc.jt^e  +  (2c.^-a£2)te  +  (£2~“£-3^® 

a^£jt^e  -  (4a£j-a^£2) te  ^  (2£j~2a£2  +  a^£2)e 


.  1,  2o,«o.o.  0,0  .  o  .  ..j 

where  £j  =  ^“^2  ^ ~  ^2'  %  ~  -1*  Associated 


histories  for  S!!»  £  are  in  Figures  3.4(a)  through  (d) 


SECONU'j 


SECONDS 


(c)  Angular  Velocity  (u) 


(d)  Applied  Torque  (]i) 


igure  3.3.  Stabilization  and  reorientation  aininizing  performance 
measure  ,  using  momentum  transfer  wheels. 
n=  1st  component,  v=  2nd  component,  0=  3rd  component 
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(a)  Attitude  Paraiceter  (y  )  (b)  Reduced  Euler  Parameters  (x) 
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(c)  Angular  Velocity  (u)  (d)  Applied  Torque  (^i) 


Figure  3.4. 


Stabilization  and  reorientation  minimizing  performance 


: 


(a)  Attitude  Parameter  (y  ) 


(b)  Reduced  Euler  Parameters  (y_) 


(c)  Angular  Velocity  (w) 


(d)  Applied  Torque  (x) 


Figure  4.2.  Detnmbling  and  orientation  minimizing  performance 
measure  with  t,  =  10  seconds. 

D=  1st  component,  ^  =  2nd  component,  0=  3rd  component 
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- — -  (1  -  ILilh  =  10-83  secs 

4  llxjll 

which  is  less  than  the  50  second  specified  final  time.  But  choosing 
i-l  “  10  secs,  the  detnmbling  and  orientation  maneuver  can  be 
accomplished  (see  Figure  4.2).  At  the  input  is  returned  to  zero 

and  the  spacecraft  maintains  the  desired  attitude. 

Using  (20)  to  determine  an  upper  limit  for  maneuver  time  results  in 
an  unnecessarily  restrictive  bound  on  t^^  The  specific  detumbling 
problem  defined  earlier  can  in  fact  be  accomplished  if  t^  chosen  to 

be  45  secs,  (see  Figure  4.3).  Finding  this  longer  acceptable  maneuver 
time  allows  torque  magnitudes  to  be  significantly  lowered. 

An  alternative  to  attempting  to  detumble  and  orient  the  spacecraft 
in  the  same  maneuver  is  to  separate  the  problem  into  the  two  distinct 
tasks;  first  detnmble  the  spacecraft,  then  reorient  to  the  desired 
attitude  in  a  rest-to-rest  task.  One  strategy  for  detumbling  is  to  set 

Pj  =  =  q2  =  =  0  and  allow  p^  in  the  performance  index 

(16),  which  yields  a  step  function  for  the  acceleration  in  the  linear 
mode  1 , 

u(t)  =  -  i,(0) . 

1  ^ 

Integrating  this  expression  twice  gives 
-2^1^  =  (1  -  ^2^0) 

.2 

ij(t)  =  Xj(0)  +  (t  -  ^) 


Applying  the  constraint. 
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Given  tbe  system's  initial  state,  and  z^*  and  the  desired 
reorientation  time  t^,  substitution  in  (19)  to  check  if  the  inequality 
is  satisfied  is  easily  accomplished.  However,  failure  to  satisfy 
inequality  (19)  does  not  mean  the  detnmbling  and  orientation  cannot  be 
accomplished  in  the  time  allotted.  In  fact,  the  implication  is  that  too 
long  a  time  has  been  specified.  This  can  be  seen  by  rearrangement  of 
the  inequality 


Now  the  bound  is  on  t^.  If  t^  as  specified  satisfies  this  relationship, 
then  the  maneuver  can  be  performed  using  the  control  associated  with 
performance  measure  If  does  not  satisfy  (20),  a  new  maneuver 
time  tj  <  tj.  is  defined  for  which  the  inequality  will  hold.  Final  time 
tj  is  replaced  in  the  optimization  solutions  by  tj  and  the  desired 
maneuver  is  performed.  From  t  =  tj  to  t  =  t£-,  zero  input  is  applied,  the 
spacecraft  attitude  does  not  change,  and  the  required  orientation  is 
achieved  in  (less  than)  the  specified  time. 

This  approach  is  demonstrated  in  solving  the  problem  introduced 
earlier  in  this  section.  Inequality  (20)  (or  (19))  is  not  satisfied 


with  the  specifications  given.  Fvalnating  the  right-hand  side  of  (20) 
gives 


(c)  Angular  Velocity  (u) 


(d)  Applied  Torque  (x) 


Figure  4.1.  Detunbling  and  orientation  minimizing  performance 
measure  with  t-  =  50  seconds. 

□  =  1st  component.  2nd  component.  o  ••  3rd  component 


becomes  0,  £  becomes  onbounded,  aod  the  simulation  program  halts. 

Investigating  equation  (18)  further,  for  4  0,  sufficient 

conditions  on  and  x^  can  be  determined  under  which  the  singularity 

will  be  avoided.  From  the  properties  of  a  norm. 


Ui(t)||  i  |-2t^^t^-t)^  +  3t"^(t^-t)^|-|UJ 


.211101 


+  I  -t“^t^-t)^  +  t/(t^-t)^|-||i^|| 


2|  I  I  0| 


As  shown  earlier,  the  polynomial  multiplying  x^  in  (18)  decreases 
monotonically  from  1  to  0  as  t  runs  from  0  to  t^.  A  few  steps  of 

calculus  determine  that  the  polynomial  multiplying  x^  is  non-negative 
over  the  time  interval,  and  has  a  maximum  of  4tj/27  at  t=tj/3.  These 
results  give 


max  (I  -  2t'^(t^-t)^  +  3t^^(t^-t)‘^| 


U;||)  =  lUJII  .  Oitlt^ 


max  (I  -t~^tj-t)^  +  t"^(t^-t)^|  •  lli^ll)  =  ^t^  lli^ll  .  0  <  t  i  t^ 


max  lUj^t) 


i  ii^ii  .  oi,  itj. 


The  last  expression  leads  to  a  bound  on  llx^ll  terms  of  ||x°j|  which 
will  guarantee  ||x^(t)||  <  1  for  the  duration  of  the  maneuver. 
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This  last  expression  is  non-positive  for  0  ^  t  <  ®  if  a  0.  As 
discussed  in  Section  3.2.4,  a  is  the  cube  root  of  a  non-negative  number 
and  is  therefore  non-negative. 

Choosing  an  inertial  frame  to  define  Xj(t^)  =  0  may  conceivably 

result  in  ||i°||  =  1  (equivalently,  =  0).  The  singularity  is  then 
excited  at  t=0.  Dividing  the  maneuver  into  two  rotations,  and  defining 
the  inertial  frame  independently  for  the  first  rotation,  allows  the 
desired  reorientation  to  be  accomplished.  With  this  accommodation  of 
the  singularity  problem  at  the  initial  time,  any  rest-to-rest  maneuver 
may  be  accomplished  without  going  to  0. 

Such  is  not  the  case  when  the  spacecraft  is  initially  tumbling. 
Consider  the  problem  of  detumbling  and  orienting  a  spacecraft  with 
external  thrusters  when  the  initial  states  are 
X(0)  =  (0.46180,  0.19167,  0.79987)^  and  ^(0)  =  (0.075,  0.075  0.075)^, 
and  the  desired  final  time  is  50  secs.  Figure  4.1  shows  the  simulation 
results  for  this  problem  using  'performance  measure  .  At  t=4.65  sec 


•  $ 
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jj<t)  .  ij 


*  4- 


(18) 


Since  u(0)  =  0  implies  =  ^2^^^  ~  r(l.(0) )«!)(0)  =  0,  the  second  term  on 


the 

right-hand 

side  is  zero 

for 

all 

0  1  t  i  t^. 

With  this 

simplification. 

x^(t)  =  p(t)i° 

where 

the 

polynomial 

p(t) 

=  -2t“^(tj-t)^ 

+  3t“^(t^-t)^ 

has 

the 

properties 

p(0)  =  1  and 

p(t^)  =  0.  Consider  the  time  derivative  of  p(t). 


=  6t'^(tj-t)^  -  6t~^(tj-t) 

=  6t^^(tj-t)  I(tj-t)  -  t^] 

=  6t~^(t^-t)  [-t]. 

The  three  terms  outside  the  bracket  are  all  non-negative  and  the  term 
[-t]  is  non-positive  for  0  t  ^  t^.  Therefore,  the  derivative  of  p(t) 
is  non-positive  and  p(t)  is  monotone  decreasing  for  the  duration  of  the 
maneuver.  Because  of  the  nature  of  this  polynomial, 

=  |p(t)|-  lli^ll  i  |p(o)|-  IlijM  =  IlijM  for  all 

0  i.  f  1  tf*  If  the  initial  orientation  of  the  spacecraft  does  not 

correspond  to  =  0,  or  equivalently  ||l°||  =  1,  the  singularity  will 
not  be  encountered  in  this  rest-to-rest  maneuver.  An  argument  similar 
to  that  just  presented  can  be  developed  for  the  remaining  performance 


measures . 
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This  chapter  addresses  the  problems  of  singularity  avoidance  and 
variation  of  parameters.  Control  laws  under  which  ^  ^ 

0  ^  t  ^  t^  are  developed.  Compensation  for  variations  from  the  modeled 
inertial  properties  is  introduced.  Three  approaches  to  the  control  of 
spacecraft  attitude  with  torque  actuator  constraints  are  discussed  in 
Chapter  Five.  Simulations  demonstrate  the  capability  to  deal  with  these 
three  difficulties. 

4.1.  Singularity  Avoidance 

In  determining  conditions  under  which  the  transformation 
singularity  will  not  be  excited,  benefits  are  realized  from  dividing  the 
spacecraft  attitude  control  proble-n  into  two  sub-tasks.  A  rest-to-rest 
maneuver  is  characterized  by  w(0)  =  0  and  w(t^)  =  0.  The  inertial 

frame  is  chosen  as  discussed  in  Chapter  Three  to  define  the  desired 

position  parameter  as  =  0.  Oetumbling  is  characterized  by 

<o(0)  4  0,.  but  <!!,(t^)  =  0.  Initial  and  final  values  of  the  Euler 

parameters  are  determined  by  choice  of  the  inertial  frame  and  the 

control  strategy. 

The  advantage  of  this  division  of  tasks  is  realized  in  considering 
the  rest-to-rest  maneuver.  If  I  I  I  ^  ^  llx^t^ll  11x^0)11 
all  0  i  t  i  tj,.  then  will  never  be  zero  and  the  singularity  will  be 
avoided.  Since  the  transformations  developed  in  Chapter  Two  define 
Xj  =  x»  since  closed  form  expressions  are  determined  for  Xj  in 
Chapter  Three,  the  condition  on  ||x(f)|l  readily  tested. 

(Throughout  the  following  discussion.  and  x  used 

interchangeably. ) 

When  performance  measure  (as  discussed  in  Section  3.2.1)  is 
minimized,  the  parametric  attitude  evolves  according  to 


CHAPTER  FOUR 


DEALING  WITH  DIFFICULTIES  I 

Solutions  of  the  optimization  problems  in  the  previous  chapter,  and 

the  associated  simulations,  are  developed  under  a  few  assumptions  which 

limit  practical  application  to  the  spacecraft  attitude  control  problem. 

Two  constraints  on  the  Euler  parameters  are  ignored.  The  identity 
2  1 1  I |2  2  2  2  2 

7^+  II3C.M  =  ''^o'*’"'^l'^''^2'''''^3  ~  ^  physically  satisfied  at  every  instant  of 

time,  but  no  provision  for  this  is  made  mathematically.  A  singularity 

exists  in  the  input  transformation  from  the  linear  model  to  the 

nonlinear  model  when  v  =0.  This  further  constrains  the  norm  of 

0 

1.  to  I Ixl I  ^  A  second  assumption  inherent  in  the  solution  is  that 

the  plant  is  perfectly  known.  Inertial  properties  of  the  spacecraft  and 
wheels  are  the  only  plant  parameters  to  appear  in  the  mathematical 
model.  Variation  of  actual  values  of  those  parameters  from  the  values 
assumed  in  the  model  act  as  a  disturbance  and  perturb  the  actual  state 
trajectory  from  the  nominal  determined  in  Chapter  Three.  Finally, 
stated  in  the  solutions  in  the  previous  chapter  is  the  assumption  that 
input  torques  have  no  physical  constraints,  or  the  constraints  are  so 
large  that  they  may  be  ignored.  When  inputs  commanded  by  the  solution 
of  the  unconstrained  optimization  problems  exceed  actuator  capabilities, 
a  clipped  version  of  the  desired  torque  will  be  applied,  and  the 
spacecraft  attitude  history  will  not  follow  the  intended  path. 
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ignored.  Since  =  x*  this  constraint  implies  llxiCt)!!  1  1  for  all  t. 
Further,  to  avoid  the  singularity  in  the  inverse  transformation  from  u^ 
to  ji,  the  strict  inequality  llx^it)]!  <  1  must  hold.  Satisfaction  of 
these  constraints,  and  compensating  for  possible  variations  in  system 
parameters  are  discussed  in  the  next  chapter. 
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for 


in  problem 


III. 


When 


qi=q2=o,  £2  = 


r  \Fi  +  i2  “  ^  stated.  A  similar  answer  is  found  if 
and  are  set  to  zero  in  case  IV. 

A  linear  differential  equation  for  Xj^(t)  can  be  found 
system  equations  and  the  closed  loop  control  law  are  known. 
Ill,  this  differential  equation  is 


when  the 
For  case 


il(t)  +  J2  +  q^  ij(t)  +  ^  Xj(t)  =  0 

The  state  penalties  can  be  chosen  to  give  a  critically  damped  response 
curve  for  Xj(t).  Similarly  in  case  V, 


Xj  (t) 


a, -2a  +2 
1  o 


* 


ij(t) 


a^+q^  Xj(t)  =  0. 


Now  a^,  aj,  qj^,  and  q2  can  be  picked  (a^,  aj  nonzero)  such  that  the 

response  is  critically  damped.  A  significant  difference  exists  between 

the  two  cases  although  the  trajectories  look  identical.  The  system 

described  in  case  V  is  damped;  the  system  of  case  III  is  not.  Open  loop 

system  poles  for  case  V  have  negative  real  parts,  whereas  open  loop 

system  poles  for  case  III  are  at  the  origin. 

An  important  point  is  that  the  same  solutions  yielding  U2  for  wheel 

control  could  be  used  to  generate  the  linear  acceleration  u^  for 

variable  thruster  control,  when  available. 

The  purpose  of  the  five  examples  presented  here  is  to  demonstrate 

the  feasibility  of  solving  the  nonlinear  spacecraft  attitude  control 

problem  using  the  linearizing  transformation  approach.  However,  to  this 

2  2  2  2 

point,  the  constraint  on  the  Euler  parameters,  y ^'^'y 2^ 3  ~  been 
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3.3  Discnssion 

The  point  must  be  emphasized  that  closed  loop  control  is  not 
appropriate  for  cases  I  and  II.  As  t  approaches  the  final  time,  the 
feedback  gains  increase  without  bound.  Even  though  the  states  and  z^ 
are  going  to  zero  and  the  result  is  in  fact  a  bounded  acceleration  u. 
the  gains  cannot  be  implemented.  If  feedback  control  is  desired.  then 
some  scheme  to  work  around  the  end  point  singularity  must  be  found. 

The  effect  of  penalizing  |iu(t}||  in  addition  to,  or  instead  of, 
penalizing  ||]i(t)||  is  to  allow  free  choice  for  initial  and  final  values 
of  u.  This  can  be  considered  a  first  attempt  at  "ac  jleration 
shaping."  For  finite  final  time,  the  ezpression  for  u^  changes  from 
linear  to  cubic  in  time  as  this  penalty  is  added.  A  corresponding 
change  is  noted  in  the  curves  shown  for  ^  as  desired.  In  case  I,  the 
torque  is  basically  a  monotone  decreasing  function,  whereas  in  case  HI, 
it  has  the  general  shape  of  a  distorted  third-order  polynomial.  The 
input  to  the  nonlinear  system  can  be  informally  considered  a  distorted 
version  of  the  input  to  the  linear  system. 

There  is  little  difference  in  the  mazimum  torque  required  in  each 
channel  as  the  penalty  shifts  from  ||u||  to  ||u||  with  finite  final 
time.  But  for  the  infinite  time  problem,  a  decrease  of  the  order  of  60% 
of  the  mazimum  is  shown  when  penalizing  |  |u| |  •  (This  appears  to  be 
contradictory  to  Murphy's  law.  A  smoother  torque  profile  is  generated, 
the  mazimum  required  torque  is  considerably  smaller,  and  the  state 
trajectories  are  still  acceptable.) 

Earlier  the  point  was  made  that  for  infinite  time  problems,  if  the 
terminal  state  penalities  and  the  state  trajectory  penalties  are  both 
zero,  the  optimal  solution  is  n=0.  This  is  affirmed  in  the  ezpression 
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3.2.5  Case  V 

System  damping  is  introduced  in  the  final  case.  Final  time  is 
infinite  and  the  penalty  on  u  is  dropped.  Therefore  p^,  p^,  P^,  and  p^ 
are  zero,  and  performance  measure  becomes 


=  2  (2.^(t)C^x(t)  +  u^(t)u(t))dt 

o 


with  system  equations 


^1  =  ^2 


-2  “^0-1  “  ^1^  *  % 


where  a^  and  a^  are  nonzero.  Solution  of  the  optimization  problem  for 


general  a^,  a^,  q^,  and  q^  results  in 


u_(t)  =  (a„ 
2  o 


-  J 


+  (*2  “ 


aJ-2ao+2  ^ a^+q^  +  q^) 


n 


Recalling  that  a^  and  a^^  are  free  parameters  introduced  in  the 
transformation  to  an  equivalent  linear  system,  they  may  be  chosen  along 
with  q^  and  q^  to  give  the  identical  response  achieved  in  case  III. 

1  t 

Letting  q^  and  q^  represent  the  state  penalties  from  the  earlier  case, 

2  1'  2  1  I  ^  ' 

set  <l2  =  =  2  *^2  ~  *1  ~  2  ^  *^2^'  Analytic  expressions 

for  Xj^(t),  X2(t),  and  u(t)  and  simulations  for  the  variables  in  the 
nonlinear  model  are  then  the  same  as  in  example  III  and  are  not  repeated 
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yields 


llil(t)||  =  llii(O)  +  (t  -  ^)  *2(0)11 

2 

<  llii(0)||  +  It  -  ^1  •  11x2(0)11 

<  llii(0)||  +  j  tjii2(0)||  <1.  0  <  t  <  tj 


h  < 


11x2(0)11 


(1  ■  11x2(0)11) 


At  first  consideration,  this  seems  to  be  an  even  more  restrictive  limit 
on  tj  than  the  previous  strategy.  However,  since  two  separate  tasks  are 
defined,  the  flexibility  in  choosing  the  inertial  frame  can  be  used  to 
increase  the  bound  on  tj .  For  the  detnmbling  portion  of  the  maneuver, 
choose  the  inertial  frame  such  that  ijCO)  =  1.(0)  =  0,  and 
x^(0)  =  r(l.(0)  )S!l(0)  =  <0(0).  The  bound  on  t^  becomes  t^^  <  2/||<j)(o)||, 
and  with  the  initial  values  given,  tj^  <  15.40  sec.  If  t^  is  chosen  to 
be  15.0  secs,  detumbling  occurs  as  t  runs  from  0  to  15.0,  the  inertial 
frame  is  redefined  at  t=15.0,  and  reorientation  is  performed  from  t=15.0 
to  t=50.0. 

Under  the  continuing  assumption  of  unconstrained  inputs,  the 

problem  of  the  singularity  in  the  input  transformation  is  successfully 

dealt  with.  The  singularity  is  encountered  in  rest-to-rest  maneuvers 

only  as  an  initial  condition.  When  this  occurs,  the  task  can  be 

accomplished  by  dividing  it  into  two  reorientation  sub-tasks.  Bounds 

are  given  for  maneuver  times  for  either  tumbl ing-to-orientation  or 

strictly  detumbling  tasks  which  will  maintain  y  ¥  0. 

o 
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A  side  effect  of  avoiding  the  singularity  is  satisfaction  of  the 

constraint  on  the  Euler  parameters.  Since  | |xl I  ^  since  is 

2  2  2  2 

not  being  directly  controlled,  the  physical  constraint  2*"^ 2  ~  ^ 

determines  the  manner  in  which  y  will  evolve. 

'o 

4.2  Variation  of  System  Parameters 

Perfect  knowledge  of  the  spacecraft  and  wheel  inertial  properties 

is  assumed  in  the  previous  chapter.  However,  as  expendable  fuel  for 

external  t^irusters  is  used,  the  matrix  I  will  change  to  some  extent. 

o 

Variations  in  will  affect  not  only  the  system  dynamics,  but  also  the 
interface  between  the  linear  acceleration  command  generator  and  the 
torque  actuator.  If  compensation  for  these  variations  is  not  provided, 
the  actual  state  trajectory  will  not  end  in  the  desired  terminal 
orientation.  This  situation  is  illustrated  in  Figure  4.4(a)  and  (b). 

Using  the  double  integrator  linear  model,  the  command  generator 

sends  the  desired  linear  acceleration  u“  =  (x)“  to  the  interface.  The 
assumed  value  of  the  inertial  matrix  is  used  in  the  interface  to 

transform  from  u^  to  x^.  Applying  x^  to  the  system  dynamics  when  is 

indeed  known  results  in  the  desired  histories  when  the 

t 

inertial  properties  have  changed,  from  to  I^,  the  desired 

trajectories  are  no  longer  followed.  This  is  shown  mathematically  by 
comparing  (x)^  X- 


d  .-.d 
u  =  (x) 


1.  , ,  d.Td,  d  p,  d.-'-l,-.  d  d.  p,  d._-ld 
-  ^  ((w  )  u  )x  +  I  (x  (I  u  X  10  )  +  I  (x  )I  X 
4  o  o  o 


X  =  -  ^  (<!)^<j))x  +  r(x)(ijj)  ^  r(x)(ijj)  ^x^ 


but 


(a)  Modeled  and  actual  agree;  desired  trajectory  results. 


(b)  Modeled  and  actual  disagree;  perturbed  trajectory  results. 


(c)  Correction  in  the  linear  atodel  for  variation  from  modeled  I  . 

o 


Figure  4.4.  Effect  of  imperfect  knowledge  of  spacecraft  inertial 
properties  on  state  trajectories. 
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An  example  demonstrates  the  effect  of  the  change  of  inertial 

parameters  on  the  state  trajectories,  A  rest-to-rest  maneuver  is 

defined  with  initial  orientation  represented  by 

T 

y.(0)  =  (0,46180,  0.19167,  0.79987)  and  desired  terminal  orientation 
y.(t^)  =  0.  Final  time  is  specified  as  25  seconds.  External  thrusters 
are  to  be  used  to  accomplish  the  reorientation.  Inertial  properties  are 
assumed  to  be  the  same  as  in  Chapter  ITiree.  Performance  measure  is 

used  to  derive  the  desired  control  strategy.  The  nominal  trajectory 
simulation  is  shown  in  Figure  4.5.  When  a  20%  perturbation  in  moments 
and  products  of  inertia  is  experienced  in  the  system  dynamics,  but  not 
in  the  interface,  the  trajectories  diverge  from  the  desired  paths 
(Figure  4.6). 

Correction  for  perturbations  in  system  parameters  is  accomplished 

by  introducing  a  Au  to  the  desired  input  u‘*[25].  This  correction  is 

generated  by  applying  a  linear  regulator  to  the  equivalent  linear 

system's  desired  and  actual  state  histories  (see  Figure  4.4  (c)). 

Desired  linear  system  state  histories  are  supplied  by  the  command 

generator.  Actual  states  of  the  equivalent  linear  system  are  obtained 

by  transforming  the  measured  physical  states  x  M*  (Assumed  values 

of  I  are  used  in  this  transformation  since  the  actual  I  is  unknown.) 
o  o 

The  error  state  equations  for  regulator  design  are  given  by 

=  to  o]  ti]  ^  E] 

where  Ax  =  x  -  x^,  Au  =  u  -  u*^,  and  w  =  x^  -  u.  As  before,  *2=1*31 

I 

because  of  the  difference  between  I  in  the  transformation  and  1  in  the 

o  o 


dynamics . 
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A  first  attempt  at  finding  Au  assumes  w  to  be  unwanted  noise. 
Modifying  the  performance  measure  appropriately,  Au  is  determined  by 
minimizing 

J  =  j  Pj(Aij(t^))^  ^1^‘f)  +  2  (Ai2(tj)) 

^f 

+  7  /  (Au(t))^  (An(t))dt 
^  0 

subject  to  the  above  state  equations.  Solution  of  this  optimization 
problem  parallels  that  given  in  Appendix  B.l,  but  to  use  feedback  to 
account  for  noise  while  avoiding  infinite  gains  at  t=t^,  the  terminal 
state  penalties  p^^  and  p^  are  left  finite.  Actual  values  of  p^^  and  P2 
are  adjusted  to  achieve  best  results.  The  correction  to  the  desired 
linear  control  is  then 


As  =  -  P  12^1.1  ~  P22^^2 
where  p^^  =  (-^  PjP2(tj-t)^  +  Pj  ( t  ^-t ) ) /d  ( t ) , 

^^22  =  ^3  PlP2^‘f"‘^^  PjCt^-t)^  +  P2)/d(t) 
and  d(t)  =  PjP2(tj.-t)^  +  J  PjCt^-t)^  +  P2(tj.-t)  +  1. 

Simulated  results  for  this  approach  to  correction  for  variation  in 
system  parameters  are  shown  in  Figure  4.7.  Values  used  for  the  terminal 
state  penalties  are  p^  and  P2  both  10^.  Assuming  w  =  -  u  to  be  noise 
appears  not  to  have  significantly  deleterious  results. 

Although  the  evolutions  of  the  Euler  parameters  do  not  appear  to  be 
highly  sensitive  to  variations  in  inertial  properties  of  the  spacecraft, 
the  problem  is  not  insignificant.  When  the  input  is  returned  to  zero  at 
t=t^  and  no  correction  is  supplied,  the  angular  velocities  are  nonzero 


0(033';  !  SNoIQbd 


and  will  remain  so.  Not  only  is  the  desired  orientation  not  achieved, 
but  the  result  of  the  "rest-to-rest "  maneuver  is  a  slowly  tumbling 
spacecraft.  Correction  for  model  variations  as  described  gives  the 
exact  desired  terminal  state,  and  the  second  difficulty  is  overcome. 
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CHAPTER  FIVE 


DEALING  WITH  DIFFICULTIES  II 

A  general  statement  of  the  spacecraft  orientation  problem  may  be 
given  as  a  requirement  to  move  from  an  initial  system  state,  £.(0)  = 

f 

to  within  a  neighborhood  of  the  desired  final  state,  £.(t^)  =  £.  ,  in  a 
specified  time  1^.  In  Chapter  Three,  five  possible  formulations  of  an 
optimal  control  problem  were  presented  which  can  lead  to  a  theoretical 
solution.  However,  since  the  problems  are  formulated  in  terms  of  the 
transformed  variables,  commanded  torques  resulting  from  those  solutions 
may  exceed  the  maximum  that  can  be  generated.  Actuator  saturation  then 

ft  d 

occurs,  and  a  clipped  version  ^2  )  of  the  desired  (2  )  torque  is 
actually  applied  to  the  system.  This  situation  is  represented  in  block 
diagram  form  in  Figure  5.1. 

Application  of  other  than  the  optimal  control  strategy  results  in  a 
suboptimal  state  trajectory  which  may  not  meet  the  specified 
requirements  and  m. y  excite  the  singularity  Yp  “  0  interface 

between  u  and  o.  A  method  of  accounting  for  input  constraints  must  be 
developed.  An  associated  problem  is  determining  whether  to  feed  the 
actual  (sub-optimal)  system  state  or  a  synthesized  (optimal)  state  to 
the  0/2  interface.  The  synthesized  state  can  be  supplied  by  the  command 


generator,  if  desired. 
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Simplest  among  the  possibilities  of  dealing  with  the  difficulty  of 
actuator  saturation  is  to  lobby  for  increased  time.  For  a  finite  time 
optimization,  this  means  increasing  t^.  An  equivalent  effect  is 

achieved  in  the  infinite  time  formulation  by  reducing  the  magnitudes  of 
the  state  penalties  q^  and  q^.  If  the  time  limitation  is  a  hard 
constraint,  then  some  other  approach  must  be  tried. 

Typically,  input  constraints  can  be  expressed  as  ^  °max' 

i  =  1,2,3.  These  can  be  appended  to  the  optimization  problem. 

However,  introduction  of  inequality  constraints  on  the  control  variable 
significantly  complicates  the  solution,  whether  the  nonlinear  model  in 
terms  of  or  the  equivalent  linear  model  in  terms  of  (x2»i2'll^ 

used.  Finding  a  single  analytic  expression  for  the  input  is  usually 
impossible . 

In  solving  a  particular  problem,  constrained  and 
unconstrained  arcs  must  be  pieced  together  to  satisfy  all 
the  necessary  conditions.  At  the  junction  points  of 

constrained  and  unconstrained  arcs,  the  control .. .may  or  may 
not  be  continuous;  ...  [26] 

Not  only  is  the  analysis  of  the  constrained  problem  more  difficult,  but 
implementation  of  the  control  strategy  developed  may  also  be  impossible. 
Physical  realization  of  control  discontinuities  can  only  be 
approximated . 

A  further  complicating  factor  arises  when  using  an  equivalent 

linear  system  to  determine  the  inputs.  Bounds  on  c  in  the  nonlinear 
model  do  not  translate  to  corresponding  bounds  on  u  in  the  linear  model 
over  the  entire  trajectory.  To  illustrate  this  point,  a  general 
representation  of  the  relationships  between  o  and  u  is  of  benefit, 

u  =  h(i^)  +  =  h  (x)  +  D  (x^lg 
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2  =  “(£.)  +  M(i.)u  =  m  (i)  +  M  (i)u. 


With  X  =  T(£.),  the  following  equalities  are  developed. 


t  (x)  =  h(T(£.)) 


m  (i)  =  m(T(i) ) 


H  (x)  =  H(T(i.)) 


M  (x)  =  M(T(i.)) 


h(l)  =  -  mi))  mil) 


mil)  =  -  (E(i))  hil) 


H(i.)  =  (M(£.))' 


M(i.)  =  mil))' 


Comparing  this  representation  to  the  transformation  developed  in  the 
second  chapter  for  use  with  external  thrusters,  mil)  M{i.)  are  given 


1  T 

m(£)  =  Z —  (u  (D-4a  )I  v  -  I  wxw  -  a.I  w 
—  2y^  0  0  o -  1  0“ 


M{i.)  =  i^r  (r) 


T(t)  =  m(£.(t))  +  M(j^(t))Uj^(t)  . 


The  time  variable  is  explicitly  shown  here  to  emphasize  that  the 
dependence  of  a  on  u  (in  this  case,  x  on  u^^)  is  state  related,  and 
therefore  implicitly  a  function  of  time.  With  M  as  defined,  each 


component  of  x  i*  dependent  on  all  components  of  u ^ .  Bounds  on  x  are 

of  the  form  x .  I  x 

'  i'  max 


From  (21) 
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1  la^(i)l  -  l"ii(£.)Oiil  *  I"i2<t>ai2l  +  "i3'i>”l3l- 

If  maximums  could  be  determined  for  im^(£.)|  and  |M.^(j£^)|,  j  =  l,2,3,  then 

maximums  could  be  assigned  to  |u..|.  j=l,2,3,  with  some  freedom,  tc 

ij 

assure  Ir.l  would  never  exceed  t  .  There  are  two  drawbacks  to  this 
1  max 

approach,  though.  First  these  bounds  on  the  components  of  u^  could  be 
extremely  conservative,  for  as  ^  approaches  0,  ©(i.)  also  approaches  0 
and  M(i.)  approaches  I^.  The  implication  is  that  as  the  terminal  state 
is  approached,  use  of  the  full  capacity  of  the  external  thrusters  would 
be  limite'.  not  by  the  optimization  of  a  measure  of  performance,  but  by 
the  artificially  low  constraints  on  the  equivalent  linear  system  input. 

A  second,  and  overriding,  problem  is  that  of  finding  the  bounds 

over  the  entire  trajectory  on  m^{i.).  Components  of  x  are  bounded  by 

2  2  2  2  t  I 

the  requirement  Yo^^l^^2^^3  ~  However,  no  upper  bounds  on  |a)^ |  have 

been  introduced.  Without  such  constraints  on  w,  no  constraints  on 

|m^(^)|  can  be  assumed. 

Despite  the  complications  as  described,  the  advantages  of  using  the 
equivalent  linear  system  approach  to  the  control  problem  justify  trying 
to  find  a  way  to  deal  with  the  difficulty  of  actuator  saturation.  Three 
possibilities  are  discussed.  All  three  presume  development  of  an 
unconstrained  control  law.  Compensation  for  saturation,  if  any,  occurs 
when  the  constraint  is  encountered. 

To  provide  a  basis  for  comparison  of  the  different  approaches,  a 
rest-to-rest  maneuver  is  defined.  An  initial  orientation  corresponding 
to  single-axis  rotations  of  5n/2,  it/3,  and  ii/4  radians  in  a  3-1-3 


convention  is  assumed. 


The  objective  of  the  maneuver  is  to  move  the 
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spacecraft  from  the  corresponding  state  x(0)  >= 

(0.46180,  0.19167,  0.79987)^  and  li^(O)  =0  to  the  terminal  state 

X(t^)  =  S.,  i0(t^)  =  0.  Requirements  which  must  be  met  for  a  result  to 

be  considered  satisfactory  are  ly^Ctj))  i  0.01,  |o)^(tj)|  i.  0.005 

rad/sec,  i=l,2,3  with  t^’^SO  seconds  a  hard  constraint.  External 

thrusters  will  be  used  to  accomplish  the  task  with  actuator  constraints 

jr.l  S.  ^’0  Nm.  The  inertial  matrix  I  is  unchanged  from  the  previous 
1  o 

chapters.  Any  of  the  five  performance  measures  discussed  in  Chapter 
Three  may  be  used  with  penalties  chosen  to  allow  the  unconstrained 
problem  to  satisfy  the  terminal  conditions.  State  trajectories  and 
input  histories  for  the  unconstrained  problem  are  shown  in  Figures  5.2 
through  5.6. 

5.1  Uncompensated  Saturation 

If  the  reality  of  input  constraints  is  either  unanticipated  or 
knowingly  ignored,  a  policy  of  uncompensated  saturation  is  implicitly 
defined.  This  policy  is  described  in  block  diagram  form  in  Figure  5.1. 
Actuators  perform  at  their  limits  when  the  commanded  torque  is  greater 
in  magnitude  than  t  .  No  compensation  is  introduced  to  correct  for 

uftX 

the  subopt imal  torque.  Mathematically, 


- 


(l  ).  =  [m(i)  + 


if  I  (l**)  .lit 
if  I  I  >  t 

'  i'  max. 


for  iBl,2,3. 


The  resulting  state  trajectory  will  correspond  to  the  desired  state 
trajectory  only  if  the  desired  torques  never  exceed  the  actuator  limits. 
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When  open  loop  control  is  used  with  no  correction  for  saturation, 
the  error  will  be  cuainlative  throughout  the  task.  If  feedback  control 
is  desired,  the  infinite  time  optimization  must  be  used  to  avoid 
infinite  terminal  gains.  Figures  5.4,  5.5  and  5.6  show  that  for  the 
infinite  time  problems,  x  is  exceeded  by  the  commanded  torque  more  in 

fllftZ 

terms  of  both  time  and  magnitude  than  for  the  finite  time  scenarios. 

The  expectation  that  trajectories  generated  without  some 
compensation  for  input  constraints  will  not  satisfy  the  requirements  is 
confirmed  by  the  simulation  results  shown  in  Figures  5.7  through  5.10. 
Infinite  time  problem  formulation  leads  to  excitation  of  the  input 
transformation  singularity  when  the  equivalent  linear  system  is  a  series 
of  integrators.  However,  introduction  of  damping  results  in  such 
significant  improvement  that  the  specifications  are  nearly  met  (Figure 
5.11). 

5.2  Optimal  Decision  Strategy 

Use  of  a  pointwise  torque  optimization  scheme  is  suggested  by  Spong 
et  al.  [271  to  compensate  for  bounded  input  in  the  control  of  robot 
manipulators.  In  terms  of  the  robotic  application,  the  objective  of  the 
optimal  decision  strategy  is  to  minimize  a  weighted  norm  of  the 
difference  between  achieved  and  desired  joint  acceleration  vectors  at 
each  sample  time  along  the  trajectory.  Statement  of  the  corresponding 
objective  for  spacecraft  applications  requires  only  replacing  "joint 
acceleration  vectors"  by  "angular  acceleration  vectors." 

A  reference  trajectory  for  the  nonlinear  system  model,  is 
required  for  this  approach.  Transformation  of  the  solution  of  the 
unconstrained  problem  conveniently  provides  this  reference  trajectory. 
The  optimal  decision  strategy  is  then  formulated  at  each  sample  time. 


WClflNS  /  SECOND 


seconds 


SECONDS 


Rfl01«N5  /  SCCONO 


(a)  Attitnde  Parameter  (y  )  (b)  Reduced  Euler  Parameters  (y) 


(c)  Angular  Velocity  (]g) 


(d)  Applied  Torque  {.y) 


Figure  5.8.  Reorientation  state  and  input  histories  minimizing  J2 
uncompensated  for  constrained  torques. 
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5.10.  Reorientation  state  and  input  histories  minimizing 
uncompensated  for  constrained  torques. 
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and  bounds  determined  for  the  maneuver  time  which  will  guarantee 
satisfaction  of  the  inequality.  The  latter  approach  is  computationally 
cleaner  and  easier  to  implement. 

Excitation  of  the  singularity  =  0,  equivalent  to  | lxl I  “ 
the  inverse  transformation  does  not  occur  in  rest-to-rest  maneuvers 
except  where  i|x(0)||  =1.  In  such  an  instance,  division  of  the 
problem  into  two  reorientation  sub-tasks,  with  appropriate  choice  of 
inertial  frame  for  each  task,  circumvents  the  singularity.  Solution  of 
the  unconstrained  problem  when  the  spacecraft  is  initially  tumbling  can 
lead  to  a  violation  of  A  shorter  maneuver  time  t^  <  t^ 

can  be  found  which  will  avoid  the  violation.  Bounds  for  this  time  are 

dependent  on  |li°|l,  control  strategy  chosen.  Another 

approach  to  singularity  avoidance  is  to  divide  the  tumbling-to- 
orientation  problem  into  a  detumbling  task  followed  by  a  rest-to-rest 
maneuver.  Inertial  frames  may  again  be  judiciously  chosen  for  each 
task.  Step  inputs  to  the  linear  model  result  in  a  more  liberal  bound  on 
the  time  to  detumble,  and  allow  the  maneuver  to  be  accomplished  with 
lower  torque  magnitudes. 

Imperfect  knowledge  of  the  spacecraft  inertial  properties  leads  to 
a  perturbed  state  trajectory  and  nonzero  terminal  angular  velocity.  A 
control  which  was  designed  to  reorient  the  system  to  a  static  desired 
attitude  in  fact  leaves  the  spacecraft  tumbling  slowly.  Introduction  of 
a  linear  regulator  acting  on  the  state  error  corrects  for  this 
perturbation  and  allows  the  desired  terminal  state  to  be  achieved. 

Torque  input  constraints  are  difficult  to  deal  witn.  Attempts  to 
append  these  constraints  to  the  optimization  problem  eliminate  the 
benefits  of  using  an  equivalent  linear  system  for  control  design.  As 
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Spacecraft  attitude  control,  a  highly  nonlinear  problem  in  the  physical 
realm,  is  greatly  simplified  by  this  linearizing  transformation 
approach. 

A  number  of  immediate  advantages  over  previous  control  design 
techniques  are  realized.  Use  of  analytical  tools  for  linear  systems 
with  the  equivalent  model  allows  design  and  implementation  of  fast, 
multi-axial  large-angle  maneuvers.  Convergence  of  a  series  of  solutions 
is  not  an  issue;  exact  solutions  are  achieved.  Subsequent  maneuvers  of 
the  same  type  do  not  require  resolving  the  problem.  The  new  control  is 
given  by  substituting  in  the  original  solution  the  new  initial  time, 
initial  conditions,  and  final  time.  Exact  terminal  control  may  be 
realized  by  appropriate  choice  of  performance  criteria  in  the  linear 


space . 


For  the  spacecraft  attitude  control  problem,  the  inverse 


-1 

— fails  to  exist  on  the  subspace  defined  by 

I 1 1  =  1.  Since  the  geometrically  imposed  identity  on  the  Euler 

2  2  2  2 

parameters  forces  ~  strategies  must  be  developed 

which  result  in  state  trajectories  for  which  ||xj^(t)||  <  1.  This 
restriction  may  be  imposed  directly  as  a  nonlinear  state  constraint  in 
the  optimization  problem  used  to  derive  the  control.  Use  of  the  full 


set  of  Euler  parameters  would  not  lead  to  singularities.  However,  the 
state  constraint  would  become  which  is  intractable  as  a 
linear  problem.  Alternatively,  an  unconstrained  problem  may  be  solved 


CHAPTER  SIX 


CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDY 


6.1  Conclasions 

Results  of  the  work  described  here  can  be  applied  to  problems  other 
than  spacecraft  attitude  control.  Many  physical  processes  are  modeled 
by  a  vector  second  order  system 

ii  -  »,(i) 

ij  -  tjii)  ♦ 


T  T  T 

where  £,  =  .  For  systems  of  this  type,  a  general  vector  second 

order  equivalent  linear  model  can  be  defined. 


^1  =^2 

^2  =  "*0^1  -  “1*2  ^  a 


where  Xj  ~  ^2  ~  ^l'  ®  *o^l*  usefulness  of 

this  equivalent  linear  model  is  determined  by  the  nature  of 


^(f.) 


for  where  this  inverse  exists. 


the  physical  input  is 


recovered  from  the  equivalent  linear  input  by 


An  important  point  to  emphasize  is  that  reinitialization  will  not 
always  produce  the  desired  results.  A  finite  amount  of  time  is  required 
after  the  last  commanded  torque  saturation  to  allow  the  reinitialized 
control  strategy  to  compensate  for  that  limiting.  Even  then,  task 
requirements  can  be  set  so  stringent  that  the  constrained  problem  cannot 
be  implemented.  The  reinitialization  approach  to  dealing  with  input 
constraints  is  proposed  as  one  tool  with  which  exact  terminal  control 
may  be  realized. 
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Figure  5.18.  Reorientation  state  and  input  histories  minimizing 
with  reinitialization  after  torque  saturation. 
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Figure  5.17.  Reorientation  state  and  input  histories  minimizing  J2 
with  reinitialization  after  torque  saturation. 
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which  are  the  desired  or  limited  values  as  appropriate,  the  current 
states  are  transformed  to  linear  model  equivalents  £^(tj^), 
n(t^}.  The  command  generator  is  then  updated  (reinitialized)  using  the 
current  sample  time  t^^  as  t^  and  the  equivalent  linear  states  and  input 
as  initial  conditions.  Final  time  t^  remains  unchanged  throughout  the 
maneuver.  When  the  commanded  torque  no  longer  exceeds  the  physical 
limits  in  any  component,  updating  of  the  parameters  in  £(t)  ceases  until 
saturation  is  again  encountered.  Effectively,  the  optimization  problem 
has  been  reinitialized  at  a  time  when  the  unconstrained  problem  is 
applicable.  Figure  S.16  pictorially  presents  this  approach. 

The  reinitialization  of  the  optimization  problem  is  only  applied  in 
the  context  of  Cases  II  and  IV.  Penalizing  |  la.!  i  allows  matching  of  the 
torque  vector  at  times  t^,  and  Discontinuities  in  torque  profiles 

are  thereby  avoided.  Simulated  results  for  these  two  cases  are  shown  in 
Figures  5.17  and  5.18. 

ioi _ Discussion 

Disregard  of  torque  actuator  limits  can  have  catastrophic  results 
in  spacecraft  attitude  control  as  shown  in  the  trajectories  of  Section 
5.1.  Use  of  a  damped  second-order  equivalent  linear  system 
significantly  improves  the  uncompensated  performance,  but  fails  to 
satisfy  terminal  requirements.  Applying  an  optimal  decision  strategy  in 
a  path-following  scenario  gives  good  results  for  infinite  time  control. 
However,  the  exact  desired  terminal  state  is  not  achieved.  By 
reinitializing  the  solution  of  the  optimal  control  problem  at  times 
throughout  the  maneuver  where  commanded  torques  exceed  actuator  limits. 


exact  terminal  control  can  again  be  achieved  in  finite  time. 
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Figure  5.15.  Reorientation  state  and  input  histories  for  the  damped 
system  with  ODS  compensation  for  constrained  torques. 
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Figure  5.14.  Reorientation  state  and  input  histories  minimizing 
with  ODS  compensation  for  constrained  torques. 
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"optimum"  constrained  x,  but  tbe  code  must  be  fast  enough  to  converge 
between  samples  for  real  time  implementation. 

The  "Spong  ODS"  technique  for  dealing  with  bounded  inputs  is 
formulated  in  the  context  of  infinite  time  control  using  state  feedback. 
Figure  5.12  is  a  block  diagram  representation  of  this  approach. 
Simulations  shown  in  Figures  5.13  through  5.15  correspond  to  cases  III, 
IV,  and  V  in  Chapter  Three. 

Reinitialization  of  Optimization  Problem 

An  alternative  approach  to  compensating  for  actuator  saturation 
makes  use  of  the  principle  of  optimality.  A  control  policy  optimal  over 
an  interval  [O.t^]  is  optimal  over  all  subintervals  [29]. 

Closed  form  expressions  given  in  Chapter  Three  for  the  control  policies 
depend  on  the  parameters  initial  time  t^,  final  time  t^,  initial  states 

x^  and  x^,  and  for  Cases  II  and  IV,  initial  linear  model  acceleration 

u^.  (Expressions  for  arbitrary  initial  time  t^  are  given  in  Appendix 
B.)  At  any  time  during  a  maneuver,  the  current  time,  current  linear 
model  state,  and,  when  appropriate,  current  linear  model  acceleration 
can  be  introduced  into  the  optimal  control  expression,  and  the  resulting 
trajectory  is  optimal  in  the  same  sense  as  before. 

This  concept  is  applied  to  the  spacecraft  attitude  control  problem 
to  develop  a  constrained  control  policy.  The  unconstrained  problem  is 
solved  initially,  using  a  performance  measure  that  penalizes  |  |r|  |  • 
When  any  commanded  torque  exceeds  the  input  constraint,  that  actuator  is 
operated  at  its  limit.  At  each  sample  time  t^  during  saturated 
operation  of  one  or  more  of  the  actuators,  the  states 
are  measured.  Along  with  the  torque  vector  l.(t^),  the  components  of 
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I 
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•*  ij*  •dvT  ! * 

minimize 


(22) 


subject  to 


•’ll  ^  '..X- 


i=1.2.3 


where  the  weighting  matrix  S  is  symmetric  and  positive  definite,  but 
otherwise  freely  chosen  to  improve  system  performance.  (Unfortunately, 
no  guidelines  are  suggested  for  choosing  S.  In  the  simulation  results 
of  the  referenced  paper  [27],  the  authors  imply  the  weighting  matrix  is 
determined  by  trial  and  error.) 

Substituting  the  known  expressions  b>  =  I  ^(I  wxm)  +  I  and 

~  o  o  o 

^  =  I  ^(I  w^xtoi^)  +  I  in  (22)  gives  an  objective  function  in  terms 

“oo  o 

of  the  free  variable  At  instant  of  time,  tj^,  all  other 

-1 


a,  (tj.),  X 

rearranging  terms,  the  problem  statement  becomes 


parameters  (I  ,sd(t.),u  (t.),  r  (t.))  are  fixed.  Factoring  out  I  "  and 
O  X  ”■  X  X  o 


minimize  ^  (l+X^tj_))^I^^SI^^(x+y(tj^)) 


(23) 


subject  to  1^.1  i  T  i=l,2,3 

'  1*  max 

where  y.  =  lQj!i(tj.)x«(t^)  -  l^iji^(t^)xj!i'*(tj.)  -  x^(t^).  Since  is  a 

-T  -1 

symmetric  positive  definite  matrix,  the  product  SI^  is  also 
symmetric  positive  definite,  and  the  objective  function  (23)  is  convex. 
With  a  convex  objective  function,  satisfaction  of  the  Kuhn  Tucker 
conditions  is  sufficient  to  solve  the  minimization  problem  [28].  Any  of 
a  number  of  quadratic  programming  methods  may  be  used  to  find  the 
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vith  the  state  constraints,  an  easier  approach  is  to  solve  the 
unconstrained  problem  and  to  compensate  after  the  fact  for  commanded 
torques  in  excess  of  the  constraint.  An  optimal  decision  strategy 
minimizes  at  each  sample  time  a  quadratic  function  of  the  difference 
between  desired  and  achievable  angular  acceleration.  Improvement  is 
realized  over  trajectories  for  which  no  compensation  for  actuator 
saturation  is  implemented.  However,  if  exact  terminal  control  is 
required,  this  method  of  correction  falls  short.  A  technique  which 
reinitializes  the  optimal  control  problem  for  the  linear  model  whenever 
saturation  is  encountered  is  introduced  here.  With  this  approach,  the 
exact  terminal  state  may  be  achieved. 

In  spite  of  the  sophisticated  spacecraft  control  algorithms 
described  herein  and  available  elsewhere,  mission  specifications  can  be 
developed  which  none  of  the  techniques  will  satisfy.  Open  lines  of 
communication  between  the  mission  planner  and  the  control  engineer  are 
still  essential. 

6.2  Recommendations  for  Further  Study 

Mathematical  descriptions  of  upper  bounds  on  maneuver  times  for 
avoiding  the  input  transformation  singularity  are  extremely 
conservative.  This  is  evidenced  by  the  example  in  Chapter  Four  where 
the  bound  was  10.83  seconds  but  the  task  could  last  45  seconds  without 
violating  the  state  constraint.  More  liberal  maneuver  time  bounds  need 
to  be  defined.  A  possible  approach  may  be  to  form 

llxj(t)|i^  =  £^(t)£^(t)  and  find  the  maximum  as  a  function  of  ||x°ll' 

lix^li'  f£*  Applying  the  inequality  | 1 1  <  1  to  that  point  in  the 


trajectory  may  lead  to  better  estimates  of  allowable  mission  durations. 
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A  single  closed  form  solution  to  the  optimization  problems  which 

incorporates  the  state  (||xjll  <  D  input  (|o^|  constraints 

would  eliminate  the  requirement  for  determining  bounds  on  maneuver  time 

or  compensation  for  actuator  saturation.  The  state  constraint  can  be 
2  2  2 

expressed  as  ^xi^^l2^^13  ^  input  constraint  can  be  stated 

No  i  c.  where 


N' 


T 

0 

.0 


-1 

0 

0 


0  0  0  0' 

1-10  0 
0  0  1-1. 


and 


T 

c 


(1,1,1,1,1,1) . 

max 


A  state  and  input  constrained  optimization  problem  would  then  be 
formulated 


minimize 


subject  to  Xj  =  x^ 


Xi(0)  =  x" 


*2  =  "®o-l  ■  “l^2 


and 


*11  *  '12  *  '13  ■  *  <  " 


No  -  c  <0 


where  the  desired  performance  measure  is  specified  for  J  and  the 
appropriate  transformation  from  o  to  u  is  substituted  in  the  last 
constraint.  Whether  a  closed  form  expression  for  u  can  be  found  is 
questionable  at  this  time,  but  such  a  solution  would  be  extremely 
helpful . 

Until  the  problem  described  above  is  solved  in  closed  form, 
techniques  such  as  those  presented  in  Chapter  Five,  and  improvements 
thereof,  must  be  employed  to  compensate  for  bounded  inputs.  The  Optimal 
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Decision  Strategy  may  be  applied  daring  satnration  in  conjunction  with 
reinitialization  of  the  closed  form  unconstrained  solutions.  Some 
improvement  in  performance  daring  actuator  saturation  may  occur,  but  the 
cost  of  this  improvement  is  that  of  running  the  quadratic  programming 
problem  between  sample  times. 

This  cost  might  be  significantly  reduced  by  the  still  experimental 
analog  quadratic  minimization  circuitry  described  by  Chua  and  Lin  in 
[301.  A  possible  control  system  could  then  be  implemented  as  follows. 
The  reinitialized,  closed  form,  unconstrained  solutions  would  be 
generated  by  a  microprocessor.  Linear  system  inputs  would  be  digitally 
transformed  to  voltages  proportional  to  the  desired  torques.  When  the 
commanded  torques  require  actuator  saturation,  these  voltages  would  be 
sent  to  the  analog  quadratic  minimization  circuitry  which  would 
instantaneously  output  the  optimal  correction  signals.  The  corrected 
voltages  would  then  drive  the  power  electronics  of  the  actuators.  Such 
a  hybrid  approach  could  sharply  reduce  the  on-line  computational  burden 
decried  by  Chief  Scientist  of  the  Air  Force,  Dr.  Allen  Stubberud  [31]. 

Finally,  perhaps  the  most  interesting  (and  most  pressing)  extension 
of  the  equivalent  linear  approach  to  spacecraft  attitude  control  is  the 
application  to  flexible  spacecraft.  Throughout  this  discussion  a  rigid 
body  has  been  assumed.  As  solar  panels  are  added  to  provide  a  renewable 
energy  source,  however,  this  rigid  body  assumption  becomes  obsolete. 
The  models  will  be  more  complex,  and  the  interfaces  may  become  extremely 
complicated,  but  the  advantages  of  control  design  for  a  linear  model 
warrant  an  attempt  at  applying  linearizing  transformations  to  flexible 


spacecraft  control 
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APPENDIX  A 


STATE  TRANSITION  MATRIX  AND  FEEDBACK  COEFFICIENTS  FOR  TBE 
SINGLE-AXIS  SPACECRAFT  AITITDDE  CONTROL  PROBLEM  SOLVED 
IN  NONLINEAR  SPACE 


Using  the  same  terminology  as  is  used  in  section  1.3,  the  state 
transition  matrix  can  be  expressed  in  the  form 


-Pit  -Pjt 

*1®  +  V 

-Pit  -p^t 

L®1® 


-Pft 
^1®  ^  ^2® 


-Pit  -p^t 

die  +  d^e 


The  coefficients  of  t  in  the  exponentials,  Pi  and  P2>  are  determined  by 
the  moment  of  inertia  I^  and  the  state  penalty  weights  qi,  q2  by 
factoring 


s^  +  -^  JIq  +  q2  s  +  2^  =  (s+Pi)  (S+P2).  (A.l) 

o'*'*  o'* 

Coefficients  Si  through  d2  can  be  determined  by  partial  fraction 
expansion  while  the  state  transition  matrix  is  in  the  Laplace  domain. 
All  the  above  coefficients  may  be  either  real  or  complex. 

Torque  input  for  this  problem  is  determined  by  feedback  of  terms  of 
successively  higher  order  in  the  state  variables, 

T  =  [biY  +  62“^  **’• 
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The  first  and  third  order  coefficients  are  listed  here. 


-- 


+ 


■  “  21  ^4p,  <aidi+3ajbjCj)  +  3p  +p  <®i^2'^^“l®2‘^l'^^®l*’l®2 

O  1  12 


+  SaJb^Cj+dajSjbjCj)  +  (3aJa2d2+3aja^dj+3aJb2C2 


+  6aja2bjC2+6aja2bjC2+6aja2b2Cj+3a2bjCj) 


^  4^  <3a2b2C2+a^d2)l 


<2aJcjd2+2aJc2dj+4aja2Cjdj 

X  12 


+  4ajbjCjC2+2a2bjCj+2a^b2cJ)  +  2pj+2p2 


+  4aja2C2dj+2a2Cidj+2a2b2cJ+482bjCjC2+4ajb2CjC2+2ajbjC2) 


+  ;-4^(2a2cjd2+2a2c2dj+4aja2C2d2+2a2bic2+2ajb2c2+4a2b2CiC2) 


4^  <2a2‘>2'*2^2a2b2C^)] 
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'4  -  -  i-h’'2  '4^  *  3?^  '3-?»,32*3.Jbjd,46.j.jbjd, 


*  3.jbJoj43.jb5o,«.^b,bjC,)  4  (3.5bjd246.j.2bidj 


+  6aja2b2dj+3a^bjdj+3a2bjC2+6ajbjb2C2+6a2bjb2Ci+3ajb 


+  p  ;3~  3a2b^d2+3a2b2dj+6a2b^b2C2+3ajb2C2 

X  2 


+  3a2b^Cj)  +  ^  (3a^b2d2+3a2b^C2)] 


2  ^  3Pj+P2  (2aJ‘li‘^2+2aja2dJ 


2  2 

+  4ajbjCjd2+4ajbjC2dj+4ajb2Cjdj+482bjCjdj+2bjCjC2+2bjb2Cj) 
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APPENDIX  B 


LINEAR-QDADRATIC  OPTIMIZATION  SOLUTIONS  FOR 
LINEAR  SYSTEM  EQUIVALENTS  OF  THE 
SPACECRAFT  ATTITUDE  CONTROL  PROBLEM 

Detailed  solutions  for  the  five  examples  presented  in  Chapter  3  are 
given  here.  The  general  problem  statement  is 

minimize  J  =  +  /  {^x'^(t)Qs,(t)  +  ^^{t  )n(t )  )dt 

o 

subject  to  x,(t)  =  Fi(t)+Ga(t),  x(0)  =  x° ,  i(tj)  =  0. 

Solution  is  by  use  of  Pontryagin's  Minimization  Principle  [32],  and 

-1  T 

gives  u  =  -  R  G  Px 

where  P  =  -  PF  -  F^P  +  PGR”^G^P  -  Q,  P(t^)  =  P^. 

Since  the  system  equations  in  each  case  are  decoupled,  the  problems  are 

T 

solved  for  the  scalar  system,  i.e.,  x  =  ,  u  =  u. 

An  initial  time  other  than  zero  will  only  affect  the  determination 
of  coefficients  in  the  analytic  expressions  for  the  state  variables.  A 

f 

desired  final  state  x  other  than  zero  can  be  achieved  by  replacing 

f 

x(tj)  by  x(tj)  -  X  in  the  finite  final  time  problems,  or  by  replacing 

£(t)  by  x(t)  -  x^  in  the  infinite  final  time  problems.  The  five 
examples  presented  in  Chapter  3  are  identified,  in  the  same  order,  in 
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App.  B.l  through  B.5  by  whether  final  time  is  finite  or  infinite, 
whether  the  control  penalty  is  on  ilp.ll  or  |  |u.|  | ,  and  whether  the  state 
equations  describe  an  integrator  or  damped  system. 
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APPENDIX  B.l 

FINITE  TIME,  PENALTY  ON  | |a | | .  INTEGRATOR  SYSTEM 
The  problem  is  stated 


...  T  1  2,^  , 

minimize  J  ~  2 

+  i  /  u^(t)dt 
‘o 

subject  to  ij  =  *2 

ij(0)  =  Xj 

Xf (tf )  =  0 

i2  =  u 

^2(0)  =  X2 

X2(tf)  =  0. 

System  and  penalty  matrices  are  defined 

-1:  ;i.  -Cl. 

,  R  =  1. 

Solution  of  the  matrix  Ricatti 

equation 

p  =  ..pF  _  F^P  + 

pgr"^g^p  -  c 

i,  p(tj)  =  p^ 

is  simplified  by  eliminating 

Q,  since 

it  is  0,  and 

finding  the 

differential  equation  for 

Z  =  (P“ 

^).  Using  the 

relationship 

2  =  -  P'^fp"^, 

2  =  FZ  +  ZF^ 

T 

-  GG^  , 

Z(t^)  =  P"^ 

Since  exact  terminal  control 

is  desired 

,  Pj  and  P2 

are  infinite. 

-1  P/Pl  °  1 

therefore  P^  =  = 

[:  :]■ 

I  11  12j 

Defining  Z  =  ,  three  differential  equations  must  be  solved, 

ri2  *22J 
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APPENDIX  B,5 

INFINITE  TIME,  PENALTY  ON  |  |a|  | ,  D.'’"^ED  SYSTEM 

The  problem  is  stated 

00 

12  2  2 
minimize  J  =  2  +  n  (t))dt 

o 

subject  to  Xj{0)  =  x° 


System  and  penalty  matrices  are  now 


1  • 

0 

It 

n 

.  P^  =  0,  Q  = 

Ql  0- 

1 

0 

1  1 

1' 

f 

^2, 

and  the  algebraic  equations  resulting  from  the  Ricatti  equation  are 

P12  =  *  ^1 

pL  =  -2»lP22-^2Pi2+q2 

Pll  ""  *1^12  ^0^22  ■*■  P12P22' 


or 


'12 


P22  “ 


Pll  “ 


\~2 

i  |aJ+2pj2+q2 
“lPl2'"“oP22-'Pl2P22* 
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Cj  =  +  201°  +  ii°)/a^. 
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*^3  “  *1* 


When  arbitrary  initial  time,  t^.  replaces  0  in  the  performance  index, 
equations  (B.4.2)  remain  unchanged,  but  the  expressions  for  the 
constants  become 


=  ^  fa^i, (t  )  -  2ax.(t  )  +  u(t  )]  e  ° 
xLlo  2  o  oj 

r  2  T 

=  l(-a  t  -a)x  (t  )  +  (2at  +l)x-(t  )  -  t  u(t  )  e 
L  olo  o2o  ooJ 

=  'T  r(o^t^+2at  +2)x, (t  )  -  (2at^+2t  )x,(t  )  +  t^u(t  )1e 
2L  o  o  lo  oo2o  ooJ 


As  in  problem  2,  solution  of  the  detumbling  equation  reduces  to  the  case 
with  Mull  penalized  by  letting  x^,  u,  and  u^^^  replace  x^,  x^,  and  u. 
Setting  q^  =  0,  the  differential  equation  for  z^(t)  then  becomes 


"l  ^  r  ^  ^3  *1  "  ^2  ‘j  =  0. 


If  q^  and  q^  are  chosen  such  that  q^  =  2  ^q^,  then 


Xj(t)  =  Cjte  +  c^e 


X2(t)  =  -  aCjte  +  (c^  -  o.c^)c 


2  ^  -at  2  »  -ot 

n(t)  =  a  Cjte  -  (2oCj-a  C2)e 


where 


Cj  =  -  (ax®  +  u°)/a 
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where 


APPENDIX  B.4 


INFINITE  TIME,  PENALTY  ON  | |u| | ,  INTEGRATOR  SYSTEM 
The  problem  is  stated 

minimize  ^  ~  2  ^  (u^^Nt))^)dt 

o 

subject  to  Xj  =  Xj(0)  =  x® 

x^  =  U  X2(0)  =  x^ 

u  =  u^^^  u(0)  =  u°. 


System  and  penalty  matrices  are  defined 


The  steady  state  Ricatti  equation  gives  six  nonlinear  coupled  algebraic 
equations 

2 

Pl3  = 


2 

**23 

■  * 

H 

2 

*’33 

■  * 

«3 

Pll 

"  ^13*’23 

P12 

=  P13P33 

P22 

"  *’23*’33 

"  ^13* 

Proper  choice  of  substitutions  leads  to  an  uncoupled  equation  in  P23 > 
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The  combination  of  the  two  conditions  implies  P22  ^  condition  1) 

with  pj]^  =  Pi2**22  Pi2  ^ 


Feedback  control  is  given  as 
n(t)  =  -  R"^G^Px(t) 


X2(t) 


(B.3.1) 


which  yields  a  differential  equation  for  X2(t) 

*1  *  *1  *  K  ‘1  -  »• 

Selection  of  the  state  penalities  q^  and  q2  determines  the  type  of 
response  exhibited  by  x^.  A  critically  damped  result  occurs  if 
q2  =  2  ^qj .  In  that  case 

Xj(t)  =  Cjte  +  C2e 

X2(t)  =  -  oCjte  +  (Cj^-ac2)e  (B.3.2) 

u(t)  =  a^Cjte  -  (2aCj-a^C2)e 


®1 


X 


o 

2 


where 
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APPENDIX  B.3 


INFINITE  TIME,  PENALTY  ON  ||u||»  INTEGRATOR  SYSTEM 


The  problem  is  stated 


minimize  J  =  2  "*■  (t))dt 

0 

subject  to  ij  =  Xj(0)  = 


i2  =  u  X2(0)  =  X2. 


System  and  penalty  matrices  are  defined 


■1:  :]•  -i:-.:!’  ■■■• 


Now  the  Ricatti  equation  becomes 


0  =  -  PF  -  F^P  +  PGR”^G^P  -  Q 


and  gives 


P12  = 


P22  “  ^^12  ^  ^2 


**11  “  P12P22' 


For  the  resulting  gain  matrix  to  be  positive  definite,  two  cond 


must  hold 


1)  >  0 


2)  P11P22  -  P12  >  0- 
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o  .  _^-4  o  j  „  ,^-3  o 
Cj  =  6t^  Xj  +  3t^  +  0.5t^  u 


i*^~4  o  _^-3  o  ,  ^~2  o 
-15t^  Xj^  -  7t^  x^  “  l.Ot^  u 


Cj  =  lOt^^Xj  +  At^x^  +  0.5tj^u°. 


When  as  arbitrary  initial  time,  t^,  is  used  in  place  of  0,  equations 
(B.2.2)  remain  the  same,  but  t^  is  replaced  by  (t^-t^)  in  the 
expressions  for  c^.  c^t  and  c^.  Detumbling  vith  arbitrary  terminal 
orientation  reduces  this  optimization  problem  to  precisely  the  form  of 
the  one  solved  in  App.  B.l  where  x^.  n,  and  u^^^  replace  x^,  and  n 

respectively.  Then 

Xi(t)  =  Cj  -  ^C2(t^-t)^  - 
.j(t)  -  OjCtj-t)’  * 

u(t)  =  -  3c2(t^-t)^  -  2cg(tj-t) 

and 


0.1.  0  .  _l_.2  o 

'l  ■  -1  *  2'f‘2  *  U'f' 


«i.~3  o  -2  o 
®2  “  “  ^^f  *2  “  *f 


®3  =  ^‘f  *2  ‘f  “  • 


Aga  in  note  the  reduction  by  one  of  the  order  of  the  polynomial 
expressions  for  x^,  X2>  and  u  when  the  specification  of  the  terminal 


state  is  dropped. 
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8**f“*’ 


31 


(tj-t) 


and 


720(t^-t)"^ 

360(t^-t)~'* 

60(t^-t)"^ 

P  =  Z~^  = 

360(tj-t)"'* 

192(tj-t)“^ 

36(tj-t)”^ 

60(tj-t)"^ 

36(tj-t)"^ 

9(t^-t)"^ 

Tlie  closed  loop  expression  for  is  then 


=  -  60(t^-t)"^i^(t)  -  36(t^-t)  -  9{t^-t)  ^u(t) 


A  differential  equation  for  z^Ct)  is  found 

a  1  3% 

*1  5(tf-t)"  x^  +  36(tj-t)  Xj  +  60(tj-t)xj  =  0 

which  is  again  an  Euler  differential  equation, 

Xj(t)  =  (tj,-t)^,  a  cubic  equation  for  r  is  found  yielding  r=5, 
r=3 .  Then 


Xj^(t)  =  Cj(tj-t)^  + 

X2(t)  =  -5Cj(tj-t)''  -  4c2(t^-t)^  -  3Cj{tj-t)^ 
u(t)=20Cj(tj-t)^  +  12c2(t^-t)^  +  eCgCtj-t) 

where  c^,  c^.  c^  are  found  by  evaluating  (B.2.2)  at  the 


(B.2.1) 


Assuming 
r=4,  and 


(B.2.2) 


boundary 


conditions 
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APPENDIX  B.2 

FINITE  TIME.  PENALTY  ON  ||n||,  INTEGRATOR  SYSTEM 
The  problem  is  stated 

minimize  J  =  ^  *  2  (n^^Nt))^dt 


subject  to 

*1  ^2 

Xj(0)  =  X® 

Xi(tf) 

=  0 

X2{0)  =  X2 

X2(tf ) 

=  0 

i  = 

n(0)  =  u° 

u(t^) 

=  0 

130 


2 

equation  becomes  a  scalar  equation  p  =  P  *  p(i£)  =  P2*  Then 

i  =  d(l/p)/dt  =  -(p/p^)  =  -1.  z(t£)  =  0.  The  solution  is 

z  =  (tj-t) 

P  =  (t£-t)  ^ 

U(t)  =  -  (t£-t)“^l2(t). 

Finding  analytic  expressions  for  ij,  and  n  in  the  same  manner  as 

earl ier. 

Xl(t)  =  *1  +  -  0.5t£^X2(t£-t)^ 

XjCt)  =  t£^X2(t£-t) 
u(t)  =  - 

An  interesting  result  is  that  when  the  final  orientation  is  specified, 
the  control  u  is  described  by  a  ramp;  but  when  the  final  orientation  is 
unspecified,  the  control  is  a  step  input. 


I 


I 
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+  6x^  =  0. 

This  Euler  differential  equation  is  solved  by  assiuning  a  solution  of  the 
form  Xj^(t)  =  (l^-t)^  and  substituting  [33],  An  algebraic  equation  in  r 
results  and  roots  r=3  and  r=2  are  found.  Therefore 


Xj(t)  =  Cj(t^-t)^  +  c^ft^-t)^. 

Differentiating. 

x^(t)  =  -3cj(t^-t)^  -  2c2(t^-t). 

Coefficients  c^  and  c^  are  determined  by  evaluating  x^(t)  and  ^2^^^ 
t=0  and  equating  those  expressions  to  the  initial  conditions.  Finally, 
the  open  loop  control  law  and  resulting  state  variable  histories  are 
given  by 


Xi(t)  =  Cj(t^-t)^  +  c^d^-t)^ 
x^Ct)  =  -3Cj(t^-t)^  -  2c2(t^-t) 
u(t)  =  6cj(t^-t)  +  20^ 


where 


-  o  -2  o 

Cl  -  -2t^  Xj  -  t^  X2 


,^“2  o  j  ^-1  o 
C2  =  3tj  xj  +  t^  x^. 


(B.1.2) 


If  a  detufflbling  maneuver  is  desired  and  the  final  spacecraft  orientation 
is  arbitrary,  p^  can  be  set  to  zero.  The  system  equation  is  simply 
x^  =  u  with  boundary  constraints  x^lO)  =  x^.  X2(tj.)  =  0,  and  the  Ricatti 
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^11  ^^12  ^  ® 

Zi2  =  Z22  ® 

Z22  =  -1 
This  results  in 

2 .  [  i<T'>'  -  2'‘f-‘>^ 

-  <tj-t) 

and 

12(tj-t)"^  6(tj-t) 

6(tf-t)~^  4(t^-t) 

Then  the  closed  loop  control  law  is 
u(t)  =  -  R"^G^P(t)i(t) 

=  -6(t^-t)“^Xj(t)  -  4(tj-t)"^X2{t) . 

From  the  system  equations, 

•  •  —5  —  1 

Xj  =  u  =  -  6(tj-t)  Xj  -  4(t^-t)  X2 

a  differential  equation  for  Xj(t)  is  found 

x'l  +  4(tf-t)“^ij  +  6(tf.-t)“^Xj  =  0 


or 


2 


Since  one  desired  result  of  introdncing  and  a^  in  the  model  is  to 
achieve  a  stable  system,  both  must  be  nonnegative.  A  matrix  Ricatti 


equation  has  one  symmetric  positive  definite  solution, 


Positive 


definite  for  P  implies  Pjj^P22  ~  ^12  ^  smaller  the  magnitude  of 
Pj2»  more  positive  the  determinant  of  P.  Therefore,  the  positive 
definite  solution  of  P  must  be  the  one  associated  with  the  smaller 


magnitude  of  V^2‘  ^12  ~  ”^o  * 


NiV*!!- 


From  the  discussion  in  App. 


B.3,  P22^®*  hence,  P22  =  +  ^®i+2pj2+q2*  these  results,  the 
expression  for  p^^  can  be  found,  but  that  expression  is  not  needed  in 
the  solution  of  this  problem.  Positive  definiteness  has  already  been 


guaranteed,  for  Pj^2  P22  suitably  chosen  (and  are  both  positive), 
and  the  expression  for  Pj^j  only  involves  nonnegative  terms.  Therefore 
P2j>0.  The  closed  loop  control  law  is  given  by 


u(t)  =  -  R'^G^Pi(t) 


(B.5.1) 


' 


ao+qi)xi(t)  +  (a^ 


which  is  independent  of 

This  control  law  yields  the  differential  equation  for  x^Ct) 


3 


